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ABSTRACT 


For  an  Energy  Management  System  (EMS)-controlled  microgrid  that  uses  a 
single-phase  voltage  source  inverter  (VSI)  configuration  to  supply  power  for  AC  loads 
from  DC  energy  storage  devices,  the  DC  link  connecting  the  VSI  to  the  DC  energy 
storage  devices  experiences  a  voltage  ripple  that  is  a  second-order  harmonic  of  the  AC 
frequency.  When  the  EMS  has  a  battery  bank  as  the  only  DC  energy  storage  device,  the 
DC  link  voltage  ripple  causes  a  second-order  hannonic  ripple  in  battery  current  that  leads 
to  a  reduction  in  battery  life.  Supercapacitors  (SCs),  which  are  capable  of  operating 
through  many  more  cycles  than  batteries,  can  be  added  to  the  EMS  to  fonn  a  battery  and 
SC  Hybrid  Energy  Storage  System  (HESS).  The  EMS  can  then  be  used  to  filter  out  the 
second-order  harmonic  ripple  in  battery  bank  current  and  provide  this  current  solely  from 
the  SC  bank,  extending  battery  life  and  reducing  system  operating  costs. 

The  design  of  a  computer  simulation  of  an  EMS  providing  power  to  a  microgrid  is 
discussed  in  this  thesis.  Simulations  provide  theoretical  results  and  laboratory  results 
compare  well  with  the  theoretical  results,  showing  that  SCs  can  be  used  to  remove  the 
second-order  hannonic  current  ripple  from  the  battery  bank. 
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I. 


INTRODUCTION 


The  Navy  Shore  Energy  Program  lists  energy  bills  as  “the  single  largest  cost  for 
Navy  installations,  reflecting  about  28%  of  Navy’s  shore  budget”  [1],  The  importance  of 
efficient  production  and  use  of  energy  throughout  the  Department  of  Defense  (DOD), 
exacerbated  by  the  increasingly  tight  fiscal  constraints  felt  throughout  all  elements  of  the 
president’s  Budget  in  recent  years,  has  led  the  DOD  to  focus  heavily  on  energy 
management.  As  a  direct  result  of  these  increased  concerns,  the  office  of  the  Deputy 
Assistant  Secretary  of  the  Navy  (Energy)  was  established  along  with  several  energy 
organizations  such  as  the  Department  of  the  Navy’s  (DON)  “Task  Force  Energy”  and  the 
United  States  Marine  Corps  (USMC)  “Expeditionary  Energy”  program.  Together,  these 
groups  strive  to  achieve  the  energy  goals  set  forth  by  the  Secretary  of  the  Navy 
(SECNAV)  and  laid  out  in  the  Office  of  the  Chief  of  Naval  Operations  (OPNAV) 
instruction  4100.5E,  “Shore  Energy  Management,”  including  the  plan  to  provide  50%  of 
all  shore-based  energy  requirements  from  alternative  sources  by  the  year  2020  [2]. 
Adding  alternative  energy  sources  and  multiple  energy  storage  methods  to  local  grids  on 
installations  will  increase  our  energy  efficiency  and  power  reliability,  but  this  will  require 
more  advanced  energy  management  techniques  to  handle  the  distributed  resources  (DR) 
providing  power  to  a  grid. 

A.  PURPOSE 

Energy  Management  System  (EMS)  technology  has  gained  much  popularity  in 
recent  years  [3],  [4],  [5],  [6],  particularly  in  allowing  the  control  of  DR  that  are  part  of  a 
microgrid.  A  microgrid  is  “a  cluster  of  micro-sources,  storage  systems  and  loads  which 
presents  itself  to  the  grid  as  a  single  entity  that  can  respond  to  central  control  signals”  [7], 
where  the  micro-sources  can  include  sources  such  as  micro-turbines  and  battery  banks. 
Among  the  many  benefits  provided  by  an  EMS  is  the  ability  to  allow  efficient  and 
reliable  sharing  of  grid  loads  by  several  disparate  power  sources,  which  allows  alternative 
energy  sources  such  as  solar  and  wind  energy  to  provide  power  to  installation  loads 
alongside  more  traditional  power  sources  such  as  diesel  generators.  The  EMS  also  allows 
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for  the  concurrent  use  of  several  different  energy  storage  methods,  maximizing  an 
installation’s  energy  efficiency  and  cost  effectiveness. 

The  use  of  traditional  batteries  as  the  sole  energy  storage  method  in  a  microgrid  is 
plagued  by  cost  and  reliability  issues.  Batteries  are  often  quite  expensive  to  purchase  and 
replace,  and  there  are  many  factors  such  as  temperature,  number  of  discharge  cycles,  and 
depth  of  discharge  that  can  negatively  affect  battery  service  life.  For  an  Energy 
Management  System  (EMS)-controlled  microgrid  that  uses  a  single-phase  voltage  source 
inverter  (VSI)  configuration  and  a  battery  bank  for  energy  storage,  the  DC  link 
connecting  the  VSI  to  the  battery  bank  experiences  a  voltage  ripple  that  causes  a 
corresponding  ripple  in  battery  bank  current,  leading  to  a  reduction  in  battery  life  and 
increased  system  operating  costs.  This  DC  link  voltage  ripple,  which  is  a  second-order 
hannonic  of  the  AC  load  frequency,  is  caused  by  the  use  of  the  VSI  to  convert  DC  to  AC 
and  will  result  in  a  second-order  harmonic  ripple  in  battery  current  unless  the  EMS  is 
configured  to  filter  out  this  ripple  current  from  the  battery  bank  and  provide  it  from 
another  source.  Supercapacitors  (SCs)  can  be  added  to  the  EMS  to  form  a  battery  and  SC 
Hybrid  Energy  Storage  System  (HESS).  The  EMS  can  then  be  used  to  filter  out  the 
second-order  harmonic  ripple  in  battery  bank  current  and  provide  this  current  solely  from 
the  SC  bank,  extending  battery  life  and  reducing  system  operating  costs.  As  stated  in  [8], 
electric  double-layer  capacitors  known  as  SCs  have  desirable  short-term  energy  storage 
and  power  properties  when  compared  with  batteries,  making  them  an  attractive  option  for 
hybrid  energy  storage  systems  (HESSs).  While  SCs  have  a  lower  energy  density  than 
batteries,  they  have  a  higher  power  density  than  batteries  and  can  provide  very  high 
output  power  (on  the  order  of  hundreds  of  kW)  for  seconds  at  a  time  [8].  In  addition  to 
their  high  power  density,  SCs  have  an  increased  life  cycle  and  can  operate  through  many 
more  charge  and  discharge  cycles  than  traditional  batteries  (more  than  100,000  [8]) 
before  requiring  replacement. 

In  this  thesis,  one  method  of  employing  an  EMS  to  provide  and  control 
installation  power  using  a  battery  and  SC  HESS  is  covered,  where  the  system  is 
configured  as  a  parallel  power  source  with  the  main  grid.  The  focus  is  on  showing  how 
the  use  of  a  battery  and  SC  HESS  reduces  battery  ripple  current,  which  can  maximize 
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battery  lifetime  and  ultimately  reduce  the  amount  of  money  that  the  DOD  spends  on 
energy  for  its  installations  each  year.  A  computer  simulation  of  the  grid  providing  power 
to  a  microgrid  via  an  EMS  and  employing  a  battery  and  SC  HESS  clearly  shows  the 
functionality  of  the  system.  The  theoretical  results  of  these  simulations  are  confirmed  by 
laboratory  experimentation. 

B.  RESEARCH  OBJECTIVES 

The  primary  research  objective  of  this  thesis  is  to  show,  using  both  simulation  and 
experiment,  that,  when  using  an  EMS  on  a  microgrid  with  a  single-phase  Voltage  Source 
Inverter  (VSI),  the  battery  bank  ripple  current  caused  by  the  second-order  hannonic  DC 
link  voltage  ripple  can  be  filtered  out  by  adding  an  SC  bank  with  a  band  pass  filter  (BPF) 
and  associated  control  circuitry  to  the  system.  In  exploring  the  primary  objective,  a 
secondary  objective  is  to  determine  if  the  battery  ripple  current  can  be  further  minimized 
by  adding  an  additional  low  pass  filter  (LPF),  leaving  essentially  only  DC  current  on  the 
battery  bus  in  steady  state. 

C.  RELATED  WORK 

There  are  several  examples  in  the  literature,  such  as  [9],  [10],  [11],  and  [12],  of  a 
battery  and  SC  HESS  used  in  grid-connected  applications  and  controlled  by  an  inverter- 
based  EMS.  While  the  benefits  achievable  by  using  this  energy  storage  scheme  are  well 
understood,  there  are  many  different  opinions  on  how  to  best  design  the  power  control 
system  to  achieve  maximum  efficiency  and  battery  life  extension.  Various  methods, 
including  the  use  of  fuzzy  logic  control  [13],  inner-loop  current  and  outer-loop  voltage 
control  [12],  and  model-predictive-based  control  [14]  have  been  covered  in  great  detail. 
The  approach  chosen  in  this  thesis  differs  from  those  found  in  literature  in  the  way  that 
the  system  model  is  built  and  analyzed  and  in  the  way  that  battery  and  SC  current  and 
power  are  controlled.  Another  focus  found  in  much  of  the  related  literature  is  to  find  a 
way  to  quantify  the  amount  of  battery  life  extension  achieved  via  the  use  of  an  EMS  and 
HESS.  In  [15],  Shen  et  al.  use  a  battery  life  degradation  model  to  optimize  battery  cycles 
in  an  electric  vehicle  with  a  battery  and  SC  HESS.  In  [16],  models  are  analyzed  by 
calculating  the  parameters  that  affect  battery  performance  and  comparing  the  differences 
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in  these  values  between  different  models  and  configurations.  The  specific  estimates  of 
battery  life  extension  afforded  by  using  the  system  that  is  described  are  not  covered  in 
detail,  but  the  observations  in  the  aforementioned  literature  and  the  general  notion  that 
reducing  battery  ripple  current  (cycles)  leads  to  an  extension  in  battery  life  are  relied 
upon. 
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II.  ENERGY  MANAGEMENT  SYSTEM  WITH 
HYBRID  ENERGY  STORAGE 


In  order  to  understand  the  operation  of  the  EMS  as  employed  in  this  thesis,  it  is 
necessary  to  learn  about  its  architecture,  including  power  converters  and  control  systems. 
The  requisite  background  information  is  provided  in  this  section. 

A.  SYSTEM  OVERVIEW  AND  DESCRIPTION 

In  this  section,  an  in-depth  overview  is  provided  of  the  EMS  that  is  used  in  both 
simulations  and  laboratory  experimentation.  The  circuit  diagram  of  a  microgrid  with  a 
HESS  and  EMS  is  provided,  and  the  simulation  circuitry  is  derived  directly  from  this 
diagram. 

1.  Physical  System  Overview 

The  EMS  used  for  this  thesis  is  a  power  electronics-based  system  designed  to 
control  and  manage  energy  sources  and  storage  on  a  microgrid.  A  microgrid  is  a  “local 
energy  grid  with  control  capability  that  can  disconnect  from  a  main  grid  (where  it  is 
connected  at  the  point  of  common  coupling)  for  autonomous  operation”  [17].  The  main 
components  of  the  EMS,  shown  in  Figure  1,  are  the  EMS  logic  circuitry,  the  primary 
controller,  and  the  Insulated  Gate  Bipolar  Transistor  (IGBT)  power  module. 

a.  EMS  Logic  Circuitry 

An  external  secondary  controller,  which  provides  an  interface  for  the  operator  to 
control  various  system  parameters  based  on  chosen  settings,  provides  input  signals  to  the 
EMS  logic  circuitry.  The  logic  circuitry  perfonns  computations  on  these  signals, 
producing  control  signals  that  are  used  to  direct  EMS  operations.  The  logic  circuitry  can 
set  the  load  current  peak-shaving  limit,  connect  or  disconnect  the  EMS  output  to  and 
from  the  grid,  perfonn  reactive  power  compensation,  shed  non-critical  loads  and 
detennine  the  direction  of  power  flow  to  manage  the  HESS.  A  complete  description  of 
the  operation  of  this  logic  circuitry  is  given  in  [18]. 
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Figure  1.  EMS  Schematic.  Adapted  from  [18]. 


b.  Primary  Controller 

The  primary  controller  module  consists  of  current  and  voltage  controllers  for  the 
different  operating  modes  of  the  EMS.  The  output  of  these  controllers  is  used  to  produce 
a  signal  that  acts  as  a  reference  value  for  the  pulse-width  modulation  (PWM)  circuitry, 
which  gates  the  inverter  H-bridge  IGBTs  so  that  the  DC  bus  can  provide  the  required  AC 
voltage  and  current  to  the  loads.  A  DC  bus  controller  is  used  to  control  the  IGBTs  in  two 
bi-directional  buck/boost  converters  (BBCs)  and  control  the  power  drawn  from  the 
battery  and  SC  modules  as  well  as  to  allow  for  battery  and  SC  charging  when  desired.  By 
employing  BBCs  between  the  battery  and  SC  banks,  we  see  that  the  system  is  arranged  as 
an  actively  controlled  HESS.  Another  possible  configuration  is  to  directly  connect  the  SC 
bank  in  parallel  with  the  battery  hank  in  what  is  known  as  a  passive  HESS.  In  this 
configuration,  battery  and  SC  bank  power  sharing  is  controlled  by  the  battery  hank 
tenninal  voltage,  along  with  battery  and  SC  internal  resistance.  Active  and  passive  HESS 
configurations  are  shown  in  Figure  2  and  Figure  3,  respectively.  As  shown  in  [19],  the 
actively  controlled  configuration  has  several  benefits,  such  as  a  much  higher  power 
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capability,  a  larger  output  voltage  range,  better  voltage  regulation,  and  smaller  system 
weight  and  volume  (among  other  benefits),  and  it  is  the  configuration  chosen  for 
this  thesis. 


Figure  2.  Active  HESS  Configuration 
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Battery 
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Figure  3.  Passive  HESS  Configuration 


c.  IGBT  Power  Module 

The  IGBT  power  module  consists  of  four  IGBTs  for  the  H-bridge  inverter  and 
two  IGBTs  for  each  of  the  two  BBCs  between  the  DC  bus  (labeled  Vbus  in  Figure  1)  and 
the  battery  and  SC  banks,  for  a  total  of  eight  IGBTs.  The  BBC  IGBTs  receive  their  gate 
signals  from  the  primary  controller,  which  modulates  them  to  control  Vbus  and  the  charge 


7 


signal  from  the  secondary  controller  (when  battery  charging  is  desired).  The  H-bridge 
IGBTs  receive  their  gate  signals  from  either  the  grid-connected  or  islanding  mode 
controller  (depending  on  mode  of  EMS  operation),  which  modulates  them  to  provide  AC 
current  to  the  loads. 

d.  Modes  of  Operation 

The  EMS  can  operate  in  either  the  grid-connected  or  islanding  mode,  based  on  an 
input  from  the  secondary  controller.  As  described  in  [18],  the  secondary  controller 
provides  four  logical  commands  to  the  EMS: 

•  Run-This  command  is,  in  effect,  an  “On-off’  switch  for  the  EMS.  When 
Run  is  low,  the  EMS  is  not  operational.  Once  Run  is  shifted  to  high,  the 
EMS  operates  to  provide  power  to  system  loads. 

•  Source  Connect-When  Run  is  high,  the  value  of  Source  Connect 
determines  whether  or  not  the  AC  grid  is  connected  to  vac,  the  AC  output 
voltage  of  the  EMS.  If  Source  Connect  is  low,  the  source  is  not  connected. 
If  Source  Connect  is  high  and  the  source  RMS  voltage  is  sufficiently  high 
(100  Vnns),  a  thyristor  is  gated  shut  to  connected  the  source  to  vac. 

•  Current  Threshold  (CT)-This  signal  sets  the  load  current  at  which  the 
EMS  uses  the  HESS  to  conduct  peak  shaving,  reducing  the  power  peaks 
supplied  by  the  source. 

•  Charge-This  signal,  when  set  to  high,  causes  the  BBC  between  the  DC 
bus  and  the  battery  bank  to  operate  in  buck  mode,  charging  the  battery 
bank  from  the  source.  In  order  for  charging  to  occur,  the  load  current  must 
be  less  than  the  CT  (indicating  that  peak  shaving  is  not  desired). 

Each  of  the  four  logical  commands  described  above  can  be  manually  entered  by  an 
operator  or  automatically  generated  based  on  system  conditions.  Together,  they  define 
the  operating  mode  of  the  EMS. 

(1)  Grid-Connected  Mode 

In  grid-connected  mode,  Run  is  set  to  high  on  the  secondary  controller,  Source 
Connect  is  set  to  high,  and  AC  grid  voltage  is  at  least  100  Vm*.  This  allows  the  HESS  to 
share  system  loads  with  the  AC  grid,  providing  power  as  necessary  during  load  increases. 
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(2)  Islanding  Mode 

In  islanding  mode,  Run  is  set  to  high  on  the  secondary  controller  but  either  Source 
Connect  is  low  or  the  AC  grid  voltage  is  not  at  least  100  Vrms,  making  it  insufficient  to 
provide  reliable  power  to  critical  loads.  This  operating  mode  is  in  effect  when  the  AC 
grid  is  lost  and  the  EMS  provides  power  from  the  HESS  to  system  loads. 

2.  Components 

The  bi-directional  BBC,  as  depicted  in  Figure  4,  operates  as  a  boost  converter 
when  the  buck  switch  is  open  and  the  boost  switch  is  modulated  or  as  a  buck  converter 
when  the  boost  switch  is  open  and  the  buck  switch  is  modulated.  This  method  of 
connecting  the  battery  and  SC  banks  to  the  DC  bus  allows  power  flow  in  either  direction 
while  the  DC  bus  and  HESS  are  operating  at  different  voltages.  The  amount  by  which  the 
DC  bus  and  the  HESS  voltages  differ  detennines  the  desired  switch  modulation  duty 
cycle. 


Dboost 


Rbuck  L  Vboost 


0 


Figure  4.  Bi-directional  BBC 


The  operation  of  the  BBC  connecting  the  battery  bank  to  the  DC  bus  is  shown  in 
Figure  5.  When  operating  as  a  boost  converter,  the  buck  switch  (labeled  Sbucki)  is  open 
while  the  boost  switch  (labeled  Sboosti)  is  modulated.  Shutting  Sb00sti  causes  the  source  to 
store  energy  in  the  inductor  by  increasing  the  inductor  current.  Once  Sboosti  is  opened,  the 
inductor  voltage  instantly  switches  polarity  and  increases  as  necessary  in  an  attempt  to 
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maintain  the  inductor  current  constant,  which  “boosts”  the  converter  output  voltage  to  the 
desired  value  and  power  is  provided  to  the  DC  bus  at  v/,us. 


When  operating  as  a  buck  converter,  the  BBC  allows  for  the  battery  or  SC  bank  to 
be  charged  by  the  DC  bus,  despite  the  fact  that  the  DC  bus  is  operating  at  a  higher 
voltage  (200  V)  than  either  bank.  In  this  configuration,  Sbucki  is  modulated  and  Sboosti  is 
kept  open,  as  shown  in  Figure  6. 


I  source 
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Shutting  Sbucki  causes  the  DC  bus  (source)  to  store  energy  in  the  inductor  by 
increasing  inductor  current  and  providing  current  to  the  batteries  and  SCs  to  charge  the 
desired  bank.  Opening  Sbucki  allows  current  flowing  to  the  selected  bank  to  dissipate, 
reducing  charging  power  flow.  By  controlling  the  duty  cycle  of  Sbucki,  we  maintain  the 
desired  charge  rate  at  the  desired  voltage. 

a.  Supercapacitors 

Electric  double-layer  capacitors,  also  called  supercapacitors  (SCs),  have  become 
very  popular  in  recent  years.  To  understand  the  benefits  of  employing  SCs  in  an  EMS 
HESS  instead  of  conventional  capacitors,  it  is  important  to  understand  the  differences 
between  these  two  devices. 

(1)  Traditional  (Electrolytic)  Capacitor  Theory 

An  electrolytic  capacitor  is  a  device  constructed  of  a  dielectric  material  (which 
does  not  conduct  electricity  but  can  store  energy  by  becoming  polarized)  sandwiched 
between  two  electrodes.  As  a  voltage  is  applied  across  the  capacitor,  the  charge  Q  stored 
in  the  electrostatic  field  between  the  positive  and  negative  electrodes  is 

Q  =  cv,  (l) 

where  C  is  the  capacitance  and  V  is  the  voltage  across  the  SC.  The  capacitance  is 
measured  in  farads  and  given  by 

C  =  So4,  (2) 

a 

where  s0  is  the  permittivity  of  free  space,  A  is  the  area  of  the  electrode  in  m  ,  and  d  is  the 
distance  between  the  two  electrodes  in  m.  It  can  be  seen  that  the  amount  of  charge  a 
traditional  electrostatic  capacitor  can  store  is  dictated  by  the  area  of  the  electrodes  and  the 
distance  between  them.  The  energy  E  that  is  stored  in  a  capacitor,  directly  related  to  the 
amount  of  charge  that  the  device  can  store  and  the  type  of  dielectric  material  chosen,  is 

E  =  —CV2  .  (3) 

2 
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(2)  Supercapacitor  Theory 

While  SCs,  also  known  as  ultracapacitors,  are  electrochemical  capacitors  that 
store  charge  using  a  different  method  than  that  used  by  conventional  capacitors,  they  are 
governed  by  the  same  basic  charge,  capacitance,  and  energy  equations.  As  described  in 
[8]  and  shown  in  Figure  7,  supercapacitors  are  constructed  of  two  porous  electrodes, 
separated  by  an  ion-conductive  membrane,  in  electrolyte  solution. 


Figure  7.  Supercapacitor  Construction.  Source:  [8]. 


A  complex  electrochemical  reaction  occurs  at  the  boundary  between  the  pores  in 
the  electrodes  and  the  electrolyte  (where  energy  is  stored),  which  displaces  ions  that 
diffuse  through  the  electrolyte  and  the  separating  membrane  between  the  electrodes.  Due 
to  the  large  surface  area  of  the  electrodes,  SCs  can  store  much  more  charge  and  attain 
much  higher  capacitance  than  traditional  capacitors. 

There  are  many  possible  mathematical  models  used  for  an  SC,  ranging  from  very 
complex  models,  such  as  the  one  described  in  [20],  to  simpler  models,  such  as  the 
dynamic  model  described  and  used  in  [21].  The  dynamic  model  of  a  supercapacitor  used 
for  simulations  in  this  thesis,  which  includes  an  ideal  capacitor  and  an  equivalent  series 
resistance  (Rsc),  is  shown  in  Figure  8.  This  model  ignores  the  effect  of  temperature, 
frequency,  and  voltage  changes  on  Rsc  over  time,  but  since  the  SC  equivalent  series 
resistance  is  unaffected  by  the  discharge  depth,  as  described  in  [22],  this  model  is 
sufficiently  accurate  for  the  descriptions  in  this  thesis. 
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G  Rsc 

Figure  8.  SC  Dynamic  Model 


b.  Single-Phase  Voltage  Source  Inverter  (H-Bridge) 

The  single-phase,  H-bridge  VSI  has  properties  that  are  unique  to  single-phase 
systems  that  are  directly  responsible  for  the  generation  of  battery  ripple  current.  A 
description  of  this  system  is  necessary  to  understand  the  cause  and  solution  of  this 
problem. 

(1)  General  Description 

Figure  9  is  a  depiction  of  the  H-bridge  inverter  used  in  the  EMS.  The  four  IGBTs 
of  the  circuit  are  gated  by  the  pulse-width  modulated  output  from  the  grid-connected  or 
islanding  mode  controller  (depending  on  EMS  operating  mode)  to  control  the  direction  of 
power  flow  between  the  AC  grid  and  the  HESS. 


Figure  9.  H-Bridge  Inverter  Schematic 
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(2)  Second-Order  Harmonic  DC  Link  Voltage  Ripple 

A  common  problem  associated  with  the  use  of  single -phase  systems  employing 
PWM  H-bridge  inverters  is  the  introduction  of  a  second-order  hannonic  voltage  ripple  (at 
twice  the  AC  load  frequency)  on  the  DC  link  In  the  case  of  an  EMS  employing  a 
battery-only  energy  storage  system,  this  ripple  is  passed  directly  to  the  battery  bank 
current,  which  leads  to  rapid  battery  cycling  and  shortened  battery  life.  The  source  of  this 
hannonic  ripple  is  found  by  examining  the  EMS  power  output 

Pout=Vaciems  >  (4) 

where  both  vac  and  iems  are  60  Hz  AC  sine  waves  (voltage  and  current)  for  the  60  Hz 
loads.  Multiplying  these  two  60  Hz  sine  waves  produces  a  120-Hz  sine  wave  component 
of  output  power,  as  we  can  see  from  the  well-known  trigonometric  identity 


sin26 


\-cos(2Q) 

2 


(5) 


The  cos(26)  second-order  hannonic  power  term  is  reflected  across  the  H-bridge  inverter 
and  causes  a  120-Hz  ripple  in  vbus.  The  power  provided  by  the  battery  bank  is 

^batt  ^ hatrhxm  , 

where  v*a«  and  ibatt  are  battery  voltage  and  battery  current,  respectively.  The  battery  bank 
maintains  essentially  a  constant  voltage  in  the  short-term  and  is  assumed  to  be  constant 
for  the  purpose  of  simulations  in  this  thesis.  As  such,  the  120-Hz  second-order  hannonic 
ripple  is  translated  directly  to  the  battery  bus  current.  This  harmonic  increases  the  current 
seen  by  energy  storage  devices,  which  leads  to  heating  of  components  and  a  reduction  in 
battery  life.  Often,  as  described  in  [23]  and  [24],  complex  and  expensive  ripple  reduction 
circuitry  or  parallel  capacitances  are  added  to  the  DC  link  to  remove  this  ripple.  Another 
method  for  removing  this  second-order  (120-Hz)  harmonic,  which  takes  advantage  of  the 
SCs  in  the  HESS  and  the  ability  of  the  EMS  to  schedule  individual  source  currents,  is 
presented  in  this  thesis. 


B.  HYBRID  ENERGY  STORAGE  SYSTEM 

In  this  thesis,  we  focus  on  a  HESS  that  makes  use  of  lead-acid  batteries  and  SCs. 
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1. 


General  Description 


The  application  of  an  EMS  provides  the  ability  to  use  multiple  energy  storage 
methods  on  a  microgrid  by  allowing  control  of  when  and  how  power  is  provided  from 
each  storage  device.  Devices  are  selected  based  on  their  characteristics,  providing  an 
overall  system  that  best  suits  the  expected  load  profile.  Two  of  the  primary  factors  to 
consider  when  selecting  an  energy  storage  scheme  are  the  power  density  and  energy 
density  of  the  chosen  devices. 

2.  Benefits  of  Using  Supercapacitors 

As  shown  in  Figure  10,  SCs  are  capable  of  offering  a  higher  power  density  than 
batteries,  allowing  them  to  provide  power  for  large  or  fast  load  changes.  This,  in 
combination  with  the  ability  to  operate  with  high  efficiency  through  a  large  number  of 
operating  cycles  (on  the  order  of  100,000  [8]),  makes  SCs  ideal  for  use  in  a  HESS  with 
peak  shaving  and  high  pulse-power  applications. 
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Figure  10.  Simplified  Ragone  Chart  Comparing  Batteries,  SCs  and  Capacitors. 

Adapted  from  [8]. 
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SCs  have  a  higher  energy  density  than  traditional  capacitors,  so  they  are  able  to 
provide  power  for  a  longer  period  of  time  before  requiring  a  recharge.  This  allows  the 
SCs  to  shave  longer  power  peaks  than  traditional  capacitors,  again  reducing  the  cycles 
seen  by  the  HESS  batteries.  Overall,  the  combination  of  batteries  and  SCs  results  in  a 
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HESS  that  has  a  higher  overall  power  density  than  a  battery-only  energy  storage  system 
and  a  higher  overall  energy  density  than  a  battery  and  capacitor  HESS. 

a.  Peak  Shaving 

With  an  EMS  and  a  battery-only  energy  storage  scheme,  the  power  distribution 
system  can  be  configured  to  allow  the  batteries  to  supply  power  during  large  load 
changes  that  result  in  power  peaks  [18].  This  is  very  beneficial  for  grid-connected 
systems,  allowing  a  reduction  in  grid  power  consumption  during  peak  hours  (resulting  in 
cost  savings).  The  downside  of  this  design,  as  shown  in  [25]  and  depicted  in  Figure  1 1,  is 
that  the  capacity  of  commonly  used  batteries  (such  as  lead  acid  and  LiFeP04)  degrades 
over  time  as  a  function  of  the  number  and  depth  of  charge  and  discharge  cycles  the 
battery  experiences.  As  peak  shaving  occurs,  battery  capacity  is  steadily  degraded  and 
service  life  reduced  until  battery  replacement  is  required. 


Number  of  Cycles  (Times) 


Figure  1 1 .  Battery  Capacity  Degradation  as  a  Function  of  Discharge  Depth 

and  Number  of  Cycles.  Source:  [25]. 

By  adding  SC  energy  storage  to  form  a  HESS,  the  EMS  can  choose  how  much 
energy  comes  from  the  charged  SC  bank  during  load  power  peaks,  removing  these  power 
cycles  from  the  battery  and  extending  battery  life.  This  is  done  by  using  an  LPF  to 

separate  the  low  frequency  load  changes  from  the  high  frequency  load  changes  (power 
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peaks)  and  using  these  signals  to  schedule  the  amount  of  power  coming  from  each  HESS 
energy  storage  bank  at  any  given  time.  A  detailed  description  of  the  operation  of  the  LPF 
is  provided  in  Chapter  III. A. 5. 

b.  Removal  of  120-Hz  Ripple  from  Battery  Bank  Current 

In  addition  to  providing  power  during  load  power  peaks,  SCs  can  also  be  used  to 
remove  the  second-order  hannonic  current  ripple  from  the  battery  bank.  The  EMS  makes 
use  of  a  BPF  to  isolate  the  second-order  hannonic  ripple  in  DC  link  voltage  and  uses  this 
signal  to  boost  SC  current  output  so  that  the  SC  bank  acts  as  the  source  for  this  current 
ripple,  removing  it  from  the  battery  bank.  This  further  extends  battery  life. 
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III.  MODELING  AND  SIMULATION 


Modeling  of  the  EMS  system  was  completed  using  MATLAB  and  SIMULINK 
software.  The  individual  systems  that  comprised  this  model,  and  how  these  systems  were 
derived,  are  covered  in  this  chapter. 

A.  EMS  MODEL 

There  are  many  possible  toolsets  for  simulating  the  EMS  in  SIMULINK.  For  the 
purpose  of  this  thesis,  it  was  determined  that  the  simulation  results  could  best  be  obtained 
by  modeling  the  circuit  using  physics-based  equations  rather  than  by  building  the 
simulation  with  components  from  the  SimPowerSystems  block  set,  allowing  for  much 
shorter  simulation  times  without  sacrificing  the  accuracy  of  the  results.  This  method  of 
modeling  the  EMS  ignores  switching  losses,  switch  tum-on  and  turn-off  times,  and 
assumes  gate  forward  bias  voltage  =  0.0  V  and  that  the  system  IGBTs  are  always  gated 
with  the  duty  cycle  needed  to  achieve  the  desired  current  flow. 

1.  Top-Level  System 

The  top-level  system  for  the  EMS  simulation,  shown  in  Figure  12,  takes  external 
inputs  from  the  secondary  controller  and  provides  the  necessary  signals  to  the  subsystems 
to  control  EMS  operation.  In  order  to  tailor  this  simulation  to  the  focus  of  this  thesis  (the 
effect  of  adding  the  SC  bank  with  a  LPF  and  a  BPF),  only  the  grid  source  connection 
status  (Source  ON  and  Source  OFF)  is  provided  as  a  secondary  controller  external  input 
to  the  system. 

The  top-level  system  includes  an  option  to  “Turn  on”  the  EMS  by  controlling  the 
start  time  for  current  injection  into  the  AC  system  when  the  grid  source  is  connected  and 
the  EMS  is  operating  as  a  constant  current  source.  The  major  subsystems  of  this  top-level 
system  include  a  PWM  subsystem  for  the  H-bridge  inverter,  a  Proportional  Integral  (PI) 
controller  that  takes  several  system  inputs  and  produces  a  reference  DC  bus  voltage  vref, 
and  the  subsystem  that  models  the  physical  system  itself.  This  physical  system  includes  a 
battery  and  SC  HESS  connected  to  the  system  loads  via  BBCs  and  the  H-bridge  inverter. 
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Figure  12.  EMS  Simulation  Top-Level  System 

2.  PWM  System 

The  PWM  subsystem,  shown  in  Figure  13,  uses  a  reference  voltage  vref  from  the 
output  of  the  PI  controller  and  a  sine  wave  with  the  chosen  H-bridge  switching  frequency 
to  create  two  signals,  SA  and  SB.  These  signals  are  provided  to  the  physical  system, 
where  SA  and  negative  SB  are  summed  to  provide  a  pulse-width  modulated  signal  at  the 
switching  frequency  with  a  peak-to-peak  amplitude  of  two. 


Gain, 
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Figure  13.  EMS  Simulation  PWM  Subsystem 
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3.  Physical  System 

The  physical  system  simulates  the  operation  of  the  portion  of  the  EMS  shown  in 
Figure  14,  which  includes  everything  between  the  DC  bus  and  the  AC  grid  source. 
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Figure  14.  Physical  System  Schematic 


Analysis  of  the  physical  system  began  with  perfonning  two  Kirchhoff  s  Voltage 
Law  (KVL)  calculations  (KVL  1  and  KVL  2)  and  a  Laplace  domain  Kirchhoff  s  Current 
Law  (KCL)  calculation  at  Node  A  to  derive  three  equations.  The  first  equation,  obtained 
from  KVL  1 ,  is 


L  filter  R  filter 


(V) 


where  vLfii,er  is  the  voltage  drop  due  to  the  EMS  output  fdter  inductance,  vRfiher  is  the 
voltage  drop  across  the  output  fdter  resistance,  and  the  other  two  terms  ( vac  and  vems)  are 
as  shown  in  Figure  14.  The  second  equation,  obtained  from  KVL  2,  is 

(8) 


where  vLin  is  the  voltage  drop  across  the  input  inductance  and  vRm  is  the  voltage  drop 
across  the  input  resistance.  The  third  equation,  obtained  from  perfonning  KCL  at  Node 
A,  is 


C. 


ems  l ems  ^ source  ^ load 


(9) 
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where  C,„  is  the  input  capacitance  seen  by  the  source,  5  is  the  Laplace  domain  derivative 
operator,  and  the  two  other  terms  (iems  and  ii,mj)  are  as  shown  in  Figure  14.  These  three 
equations  were  used  to  build  the  SIMULINK  model  equivalent  for  Figure  14,  as  shown  in 
Figure  15. 


Figure  15.  EMS  Simulation  Physical  System 


The  pulse-width  modulated  AC  voltage  output  of  the  H-bridge  inverter  is 
simulated  directly  without  the  need  to  generate  H-bridge  gate  signals  by  subtracting  SB 
from  SA  (both  were  retrieved  from  the  PWM  subsystem)  and  performing  a  multiplication 
with  the  DC  bus  voltage  vhus  to  produce  vac.  Using  Switch  1,  the  simulation  also  affords 
the  option  of  multiplying  vref  (an  output  of  the  PI  controller)  and  (SA  minus  SB)  directly, 
which  has  the  effect  of  producing  vac  without  simulating  the  pulse-width  modulation  of 
the  H-bridge  inverter.  This  method  of  producing  vac  does  not  affect  the  simulation  output 
but  allows  a  larger  simulation  time-step  size  and,  therefore,  shorter  overall  simulation 

times.  This  method  was  used  extensively  for  this  thesis. 
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4.  HESS  and  DC  Bus  System 

This  subsystem  simulates  the  portions  of  the  EMS  shown  in  Figure  16,  which 
includes  the  Battery  and  SC  HESS  and  the  BBCs  that  connect  them  to  the  DC  bus.  In 
order  to  maximize  simulation  speed,  the  individual  switch  modulation  is  not  modeled. 
Instead,  it  is  assumed  that  modulation  takes  place  as  necessary  to  achieve  the  desired 
boost  currents  (koost  and  iscboosi),  which  are  set  by  the  boost  controller. 


i-emsdc 


Figure  16.  HESS  and  DC  Bus  System  Schematic 

The  HESS  SIMULINK  model  takes  the  EMS  output  current  iems  and  voltage  vac 
and  employs  a  boost  controller  to  derive  the  HESS  boost  currents  necessary  to  maintain 
DC  bus  voltage  VbUS- 

a.  Derivation  of  SC  Boost  Current  (iscboost)  Model 

To  create  the  model  for  SC  boost  current  iscboost.  Figure  16  was  analyzed  to  derive 
the  equations  relating  the  desired  HESS  system  current  and  voltage.  The  SC  circuit 
schematic  used  to  obtain  the  equations  for  iscboost  is  shown  in  Figure  17.  The  boost 
controller  determines  iscboost,  and  the  DC  Bus  Controller  schedules  the  SC  bank  to  provide 
iscboost  to  the  DC  bus  by  adjusting  SC  BBC  IGBT  gating. 


23 


t-scbus 


Figure  17.  SC  Converter  Schematic 


Using  the  power  conservation  law,  P;n  =  Pout,  we  see  that  the  power  flow  into 
Node  A  in  Figure  17  must  equal  power  out  of  the  same  node.  Ignoring  the  small  voltage 
drop  across  the  SC  ESR  (8.1  mfi)  for  the  purposes  of  the  simulation,  we  get 

^sdscboost  ^bus^scbus  5  (1H) 

where  all  terms  are  as  shown  in  Figure  17.  This  equation  can  be  rearranged  to  give  the 
equation  for  SC  boost  current, 


^ 'bus  l scbus 


scboost 


(11) 


b.  Derivation  of  Battery  Boost  Current  (hoost)  Model 

The  battery  bank  circuit  schematic  used  to  obtain  the  equations  for  iboost  is  shown 
in  Figure  18.  This  current  is  determined  by  the  boost  controller  and  scheduled  by  the  DC 
bus  controller  to  be  provided  to  the  DC  bus  by  the  battery  bank. 


tbbus 


Figure  18.  Battery  Bank  Converter  Schematic 
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The  derivation  of  itoost  involves  the  application  of  the  power  conservation  law  at 
Node  A,  KCL  in  the  Laplace  domain  at  Node  B,  and  application  of  KVL  as  shown  in 
Figure  18.  The  first  equation,  produced  by  a  power  conservation  analysis  at  Node  A,  is 


^ cbatthoost  ^ bushbus  5 


(12) 


where  all  terms  are  as  shown  in  Figure  18.  The  second  equation,  produced  by  performing 
KCL  analysis  at  Node  B,  is 

^ batt  ^ cbatt  ^ boost  ?  (13) 

where  all  terms  are  as  shown  in  Figure  18.  The  final  equation,  produced  by  performing  a 
KVL,  is 

~Vhatt  +  Rbatthatt  +  LS  (hat,  )  +  V cbatt  =  0  »  (  14) 


where  Rbatt  is  the  battery  internal  resistance,  s  is  the  Laplace  domain  derivative  operator, 
and  all  other  terms  are  as  shown  in  Figure  18. 


c.  Derivation  of  DC  Bus  Current 

In  order  to  connect  the  simulation  parameters  on  the  HESS-side  of  the  H-bridge 
inverter  with  those  on  the  AC  grid  side,  Figure  19  was  analyzed  to  derive  an  equation  for 
DC  Bus  current  i debus- 
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The  derivation  of  an  equation  for  i debus  involves  the  application  of  the  power 
conservation  law  at  Node  A  and  the  application  of  KCL  in  the  Laplace  domain  at  Node 
B.  The  power  conservation  law  at  Node  A  gives 


v  1  =  v,  i  ,  , 

ac  ems  bus  emsac  ‘ 


where  all  terms  are  as  shown  in  Figure  19.  The  KCL  analysis  at  Node  B  gives 


i emsdc  h:>bus  debits  ^ bus  bus  1 


(15) 


(16) 


where  Cims  is  DC  bus  capacitance  and  all  other  terms  are  as  shown  in  Figure  19.  Next, 
equations  (15)  and  (16)  were  built  into  the  SIMULINK  model  given  in  Figure  20. 


Figure  20.  EMS  Simulation  HESS  and  DC  Bus  System 


5.  Configurable  Boost  Controller  System 

In  order  to  extend  battery  bank  life,  it  is  desirable  to  provide  current  from  the  SC 
bank  for  all  but  the  lowest  (essentially  DC)  frequency  Vbus  changes  and  provide  additional 
current  from  the  SC  bank  to  remove  the  second-order  hannonic  (120-Hz)  current  ripple 
from  the  battery  bank.  The  SIMULINK  model  of  a  configurable  boost  controller  that 
receives  a  DC  bus  voltage  error  (the  difference  between  vbus  and  the  200  V  vbus_ref), 
sends  it  through  a  PI  controller,  and  determines  the  battery  bus  boost  current  iboosi 
necessary  to  drive  voltage  error  toward  zero  to  maintain  Vbus  at  200  V  in  steady  state  is 
shown  in  Figure  21.  To  allow  the  introduction  of  a  BPF  and  LPF  into  the  system  during 
simulations,  the  boost  controller  was  designed  to  allow  easy  switching  between  three 
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different  operating  modes  by  adjusting  the  value  of  two  constants  in  the  model  MATLAB 
script  provided  in  Appendix  A,  “LPFconnectionstatus”  and  “SC  connection  status.” 
The  three  operating  modes  are: 

1 .  SC  bank  disconnected 

2.  SC  bank  connected  with  120-Hz  BPF  only 

3.  SC  bank  connected  with  120-Hz  BPF  and  LPF. 


Disconnect2 


Figure  21.  EMS  Simulation  Configurable  Boost  Controller 

The  boost  controller  LPF  Laplace  domain  transfer  function  T(s),  which  defines 
the  relationship  between  the  filter  output  voltage  and  the  filter  input  voltage  as  a  function 
of  the  input  signal  frequency,  is  given  by  [26] 

T(s)  =  — — ^  ,  (17) 

S  +  (DC 

where  5  =jco  (the  square-root  of  negative  one  times  frequency  in  radians  per  second)  and 
<x>c  is  the  filter  cutoff  frequency  in  radians  per  second.  For  all  simulations  that  employ  the 
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LPF,  coc=  0.1  rad/s,  ensuring  that  only  the  lowest  frequency  load  signals  are  used  to 

create  the  battery  boost  current  signal.  The  boost  controller  BPF  Laplace  domain  transfer 
function  is  given  [26] 


T(s)  = 


Ks 


(18) 


+  5  +(Qo}0) 

where  K  is  the  filter  gain  at  mid-band,  cao  is  the  filter  center  frequency  in  radians  per 


v®oy 


second  (24071  Hz),  and  Q  is  the  filter  Q-factor  (center  frequency  divided  by  3-dB 
bandwidth),  which  defines  the  width  of  the  filter  passband.  The  simulations  in  this  thesis 
employ  an  active  BPF  with  a  gain  K  of  30,  chosen  to  be  greater  than  1.0  to  maximize  the 
removal  of  the  second  order  harmonic  ripple  from  the  battery  bank  and  a  Q-factor  Q  =  10 
to  allow  for  a  reasonably  narrow  passband. 


The  boost  converter  provides  two  output  signals  to  the  HESS  and  DC  Bus 
subsystem,  iboost  and  iscboost ■  The  DC  bus  voltage  error  is  driven  toward  zero  by  the  iboost 
current  coming  from  the  battery  bus,  and  iscboost  is  used  to  control  the  amount  of  load  that 
the  SC  bank  takes  off  of  the  battery  bank  (to  include  the  current  ripple  caused  by  120-Hz 
DC  link  ripple).  For  the  simulations  in  this  thesis,  the  iboost  and  iscboost  signals  are  driven 
directly  by  the  output  of  the  boost  controller  such  that  their  value  always  exactly  matches 
the  commanded  signals  from  the  controller.  This  allows  for  much  faster  simulations 
without  sacrificing  the  accuracy  of  results.  For  laboratory  experimentation,  by  contrast, 
these  two  boost  current  signals  are  sent  to  the  DC  bus  controller  to  modulate  the  BBC 
switches  in  order  to  achieve  the  desired  iboost  and  iscboost  values. 


Please  note  that  all  boost  controller  figures  in  this  section  depict  switches  in  the 
positions  associated  with  the  applicable  configuration,  with  wires  and  components  greyed 
out  when  they  are  not  being  used. 


a.  SC  Bank  Disconnected 

Simulations  were  performed  with  the  SC  bank  disconnected  to  show  the  worst- 
case  battery  current  ripple  that  is  present  during  EMS  operation  when  the  battery  bank  is 
the  only  available  energy  source.  In  this  operating  mode,  no  filters  are  used.  The  output 
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of  the  PI  controller  directly  sets  the  value  of  i boost,  and  iscboost  is  set  to  a  constant  0.0  A 
(effectively  removing  the  SC  bank  from  the  system).  When  configured  for  this  mode  of 
operation,  the  boost  controller  schematic  is  as  shown  in  Figure  22. 


Disconnect 


Figure  22.  EMS  Simulation  Configurable  Boost  Controller- 

SC  Bank  Disconnected 

b.  SC  Bank  Connected  with  120-Hz  BPF  Only 

To  show  the  benefits  of  adding  the  SC  bank  to  the  system,  the  boost  controller  is 
configured  to  send  the  input  voltage  error  signal  verror  through  a  120-Hz  BPF,  the  output 
of  which  is  used  to  create  the  iscboost  signal.  The  iscboosi  signal  adjusts  the  current  coming 
out  of  the  SC  bank  in  order  to  remove  the  120-Hz  ripple  current  from  the  battery  and 
place  the  load  requirement  on  the  SC  bank.  While  this  configuration  removes  much  of  the 
120-Hz  current  ripple  from  the  battery  bank,  it  allows  all  other  current  ripple  frequencies 
(such  as  those  that  are  encountered  during  a  step  change  in  load)  to  remain  on  the  battery. 
This  boost  controller  configuration  is  depicted  in  Figure  23. 
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Figure  23.  EMS  Simulation  Configurable  Boost  Controller- 
SC  Bank  Connected  with  BPF 


c.  SC  Bank  Connected  with  120-Hz  BPF  and  LPF 

While  the  120-Hz  ripple  in  battery  current  is  greatly  reduced  by  the  addition  of 
the  SC  bank  with  the  BPF,  adding  a  LPF  on  the  output  of  the  PI  controller  with  a  very 
low  cutoff  frequency  coc  is  expected  to  reduce  the  120-Hz  current  ripple  to  essentially  0.0 
A  and  allow  battery  current  to  transition  smoothly  from  its  initial  to  final  steady-state 
value  when  commanded  battery  current  changes  (such  as  during  a  load  change).  As 
shown  in  Figure  24,  where  the  boost  controller  is  configured  with  the  SC  connected  and 
both  a  BPF  and  LPF  employed,  the  signal  prior  to  the  LPF  must  be  added  to  the  output  of 
the  BPF  and  the  signal  after  the  LPF  must  be  subtracted  from  the  output  of  the  BPF  in 
order  to  ensure  that  all  of  the  error  signal  that  was  removed  from  the  iboost  signal  is  used  to 
increase  iscboost ■  This  allows  the  SC  hank  to  remove  as  much  of  the  120-Hz  ripple  from  the 
battery  bank  current  as  possible. 
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Figure  24.  EMS  Simulation  Configurable  Boost  Controller- 
SC  Bank  Connected  with  BPF  and  LPF 


B.  SIMULATION  FOR  DIFFERENT  LOAD  AND  SYSTEM 

CHARACTERISTICS 

The  simulations  performed  demonstrate  the  effect  of  using  a  battery  and  SC 
HESS  on  battery  utilization,  including  the  removal  of  the  120-Hz  ripple  from  the  battery 
bank  current.  To  make  these  effects  apparent,  three  SIMULINK  models  are  compared 
under  given  load  conditions:  a  model  that  includes  a  battery-only  energy  storage  system 
(SC  bank  disconnected),  a  model  that  includes  the  battery  and  SC  HESS  with  a  120-Hz 
BPF  only,  and  a  model  that  includes  the  battery  and  SC  HESS  with  both  a  120-Hz  BPF 
and  a  LPF.  These  scenarios  were  modeled  by  changing  the  boost  current  controller  as 
covered  in  the  Configurable  Boost  Controller,  Section  A. 5.  For  all  scenarios,  it  is  shown 
that  the  addition  of  the  SC  with  BPF  greatly  reduces  battery  ripple  current,  and  this  ripple 
can  be  driven  down  further  by  the  addition  of  a  LPF.  Additional  measurements  and 
calculations  for  both  all  simulations  are  provided  in  Appendix  B. 
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1.  Grid  Connected,  Operator  Increases  EMS  Commanded  Current 

When  the  EMS  is  in  “Grid  Connected”  mode,  with  the  grid  power  source 
connected  to  the  loads  and  EMS  circuitry,  the  EMS  operates  as  a  constant  current  source. 
The  EMS  current  is  controlled  by  settings  on  the  secondary  controller,  which  an  operator 
can  change  to  meet  desired  conditions  as  the  HESS  and  grid  source  share  the  load.  In  this 
scenario,  at  a  simulation  time  of  10.0  s,  a  system  operator  provides  input  to  the  EMS 
primary  controller  to  increase  EMS  output  current  iems  from  1.0  A  to  4.0  A.  System  load 
is  set  at  10.0  Q  for  all  three  simulations  below. 

a.  SC  Bank  Disconnected 

With  the  SC  bank  disconnected,  the  battery  bank  has  large  peak-to-peak  current 
ripples  that  increase  as  iems  is  increased,  as  shown  in  Figure  25. 


Time  (s) 


Figure  25.  Battery  Current,  EMS  Commanded  Current  Change- 

SC  Bank  Disconnected 

The  boost  current  controller  senses  DC  bus  voltage  error  ( Vbus-vrej )  and  controls 
battery  boost  current  as  required  to  maintain  Vbus  at  200  V.  The  120-Hz  ripple  in  y/„„ 
causes  a  ripple  in  the  generated  iscboost  signal,  which  transfers  the  current  ripple  to  the 
battery  bank.  Battery  peak-to-peak  AC  current  ripple  increases  dramatically  from  1.01 
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Apk-pk  before  the  commanded  current  change  to  3.95  Apk-Pk  after  the  step  change  that 
occurs  at  10.0  s. 

b.  SC  Bank  Connected  with  120-Hz  BPF  Only 

In  this  configuration,  the  boost  controller  sends  the  DC  bus  voltage  error  through 
the  120-Hz  BPF  and  uses  this  signal  to  adjust  iscboost  so  that  the  SC  bank  provides  the 
ripple  current  caused  by  the  ripple  in  Vbus  instead  of  the  battery  bank.  As  shown  in  Figure 
26,  the  addition  of  the  SC  with  a  120-Hz  BPF  greatly  reduces  the  ripple  in  battery  bank 
current  both  before  and  after  the  EMS  commanded  current  change  and  also  allows 
battery  current  to  increase  more  slowly  from  the  initial  steady-state  average  value  to  the 
final  steady-state  value. 


Figure  26.  Battery  and  SC  Current,  EMS  Commanded  Current  Change- 

SC  Bank  Connected  with  BPF 


The  battery  bank  peak-to-peak  current  ripple  is  reduced  by  approximately  89.90% 
at  all  points  in  the  simulation  when  compared  to  the  values  obtained  with  the  SC 
disconnected.  From  Figure  27,  which  is  a  close  up  view  of  the  battery  current  after  the 
commanded  EMS  output  current  change  in  Figure  26,  it  can  be  seen  that  the  120-Hz 
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ripple  of  0.40  mApk.pk  is  still  present  on  the  battery  bank  after  the  commanded  EMS 
output  current  change  despite  the  addition  of  the  SC  bank  and  BPF  in  the  boost  current 
controller. 


Time  (s) 

Figure  27.  Battery  Current  120-Hz  Ripple,  EMS  Commanded  Current  Change- 

SC  Bank  Connected  with  BPF 

c.  SC  Bank  Connected  with  120-Hz  BPF  and  LPF 

In  this  configuration,  an  EPF  is  added  to  the  boost  current  controller  to  remove  all 
but  the  lowest  frequencies  (the  simulation  uses  0.1  radians  per  second  as  the  threshold 
frequency)  from  the  battery  boost  current  signal  iboosl  .  Frequencies  above  the  threshold 
are  used  to  increase  SC  bank  boost  current  iscboost  and  remove  these  higher  frequency  load 
requirements  from  the  battery.  As  shown  in  Figure  28,  the  approximate  steady-state 
average  value  of  battery  current  before  the  commanded  current  change,  0.83  A,  does  not 
vary  significantly  from  the  value  seen  in  Figure  26  with  the  system  using  a  BPF  only. 
Similarly,  though  not  depicted  in  Figure  28  due  to  the  fact  that  ibatt  takes  much  longer  (5.0 
s)  to  reach  steady  state  in  this  configuration,  the  steady-state  average  value  of  battery 
current  after  the  commanded  current  change  is  also  the  same  as  in  Figure  26  (3.40  A). 
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Figure  28.  Battery  and  SC  Current,  EMS  Commanded  Current  Change- 

SC  Bank  Connected  with  BPF  and  LPF 

The  peak-to-peak  SC  current  increased  from  3.05  Apk_pk  to  1 1.92  Apk.pk  due  to  the 
commanded  current  change,  but  the  SC  bank  is  well-suited  to  handle  this  larger  peak-to- 
peak  current,  and  this  only  represents  an  additional  0.49  Apk_pk  increase  in  SC  current 
ripple  when  compared  with  the  “BPF  only”  scenario.  This  modest  increase  in  peak-to- 
peak  SC  current  allows  the  SC  bank  to  remove  an  additional  10.06%  of  the  120-Hz 
battery  bank  ripple  current,  bringing  the  total  removal  to  99.99%  when  compared  with 
the  “SC  Disconnected”  scenario.  A  close-up  view  of  the  steady-state  battery  current  after 
the  change  in  EMS  commanded  current  change  is  provided  in  Figure  29. 


Figure  29.  Battery  Current  120-Hz  Ripple,  EMS  Commanded  Current  Change- 

SC  Bank  Connected  BPF  and  LPF 
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The  advantage  of  adding  the  LPF  to  the  boost  controller  circuitry  is  made  clear  in 
Figure  29.  The  0.4  Apk-Pk  120-Hz  ripple  current  that  was  previously  seen  on  the  battery 
after  the  commanded  EMS  output  change  has  been  reduced  to  0.0003  Apk-Pk,  which  is 
essentially  a  complete  removal  of  this  ripple  current. 

2.  Loss  of  AC  Grid  Source  and  Transition  to  Islanding  Mode 

Under  normal  conditions,  the  EMS  operates  in  grid-connected  mode  with  the  grid 
providing  most  of  the  power  for  system  loads.  The  power  provided  to  the  loads  by  the 
EMS  is  governed  by  the  EMS  current  and  peak-power  shaving  settings  on  the  secondary 
controller,  and  the  HESS  stores  the  generated  power  from  alternate  power  sources,  such 
as  a  photovoltaic  (PV)  cell  array,  until  needed.  One  of  the  many  benefits  of  using  an  EMS 
with  energy  storage  is  the  ability  to  continue  providing  power  from  the  HESS  to  the 
system  loads  when  the  main  grid  power  source  is  lost.  With  a  battery-only  energy  storage 
system,  the  large  transients  that  result  from  a  sudden  loss  of  grid  cause  fluctuations  in 
battery  bank  current,  quickly  reducing  battery  stored  energy  and  reducing  battery  life. 
When  employing  an  SC  hank  with  a  BPF  and  LPF  as  part  of  a  HESS,  the  high  frequency 
components  of  the  system  currents  can  be  shifted  to  the  SC  to  extend  battery  life.  The 
battery  life  extension  afforded  by  the  removal  of  the  120-Hz  battery  current  ripple  will 
prove  particularly  valuable  if  grid  restoration  is  not  possible  in  the  near  future  since  the 
length  of  time  that  the  EMS  can  operate  in  islanding  mode  is  extended  and  power  for  all 
system  loads  is  provided  without  requiring  a  system  shutdown  for  replacement  of  HESS 
batteries.  In  each  of  the  three  following  simulations,  the  grid  source  is  removed  at  a 
simulation  time  of  5.0  s,  load  is  set  at  10.0  Q,  and  EMS  reference  current  is  1.0  A. 

a.  SC  Bank  Disconnected 

In  this  configuration,  the  loss  of  the  AC  grid  at  5.0  s  results  in  a  drastic  and  almost 
immediate  increase  in  battery  ripple  current  from  1.01  Apk-Pk  to  17.03  Apk-Pk  in  steady 
state,  as  shown  in  Figure  30.  The  large  battery  current  fluctuations  can  greatly  reduce 
battery  life,  and  it  is  undesirable  to  operate  on  the  batteries  in  this  way  for  long  periods 
of  time. 


36 


Figure  30.  Battery  Current,  Loss  of  AC  Grid  Source- 
SC  Bank  Disconnected 

b.  SC  Bank  Connected  with  120-Hz  BPF  Only 

As  shown  in  Figure  31,  the  addition  of  the  SC  bank  and  BPF  allows  the  SC  to 
greatly  reduce  the  battery  current  ripple  as  the  AC  source  is  lost.  Average  battery  current 
increases  from  0.83  A  to  its  new  steady-state  value  of  14.60  A  during  the  0.8  s  following 
the  loss  of  the  AC  grid.  SC  ripple  current  increases  almost  immediately  from  2.88  Apk-pk 
to  49.48  Apk-Pk  to  provide  power  for  the  increased  load  and  continue  removing  the  120-Hz 
ripple  from  the  battery  current. 


Figure  3 1 .  Battery  and  SC  current,  Loss  of  AC  Grid  Source- 
SC  Bank  Connected  with  BPF 
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In  Figure  32,  which  is  a  close-up  view  of  the  battery  current  in  Figure  31,  the  120- 
Hz  current  ripple  of  1.76  Apk-Pk  (after  loss  of  the  AC  grid  source)  can  still  be  seen.  This 
represents  an  89.65%  reduction  in  battery  current  ripple  when  compared  with  the  “SC 
Disconnected”  scenario,  which  is  essentially  the  same  reduction  obtained  by  adding  the 
SC  with  BPF  in  the  “Grid  Connected,  Operator  Increases  EMS  Commanded  Current” 
simulation. 


Figure  32.  Battery  Current  120-Hz  Ripple,  Loss  of  AC  Grid  Source- 

SC  Bank  Connected  with  BPF 

c.  SC  Bank  Connected  with  120-Hz  BPF  and  LPF 

Similar  to  the  effect  of  adding  the  LPF  in  the  “Grid  Connected,  Operator 
Increases  EMS  Commanded  Current”  simulation,  the  LPF  removes  all  but  the  lowest 
frequency  loads  from  the  battery  boost  current  signal  iboost  and  causes  all  other  load 
frequencies  to  increase  SC  bank  boost  current  iscboost,  removing  these  higher  frequency 
load  requirements  from  the  battery.  As  shown  in  Figure  33,  the  approximate  steady-state 
average  value  of  battery  current  before  the  commanded  current  change,  0.82  A,  does  not 
vary  significantly  from  the  value  seen  in  Figure  3 1  with  the  system  using  a  BPF  only. 
Similarly,  though  not  depicted  in  Figure  33  due  to  the  fact  that  it,an  takes  much  longer 
(10.0  s)  to  reach  steady  state  in  this  configuration,  the  steady-state  average  value  of 
battery  current  after  the  commanded  current  change  is  also  the  same  as  in  Figure  31 
(14.83  A). 
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Figure  33.  Battery  and  SC  Current,  Loss  of  AC  Grid  Source- 
SC  Bank  Connected  with  BPF  and  LPF 

As  shown  in  Figure  33,  the  addition  of  the  LPF  increases  SC  current  ripple  from 
3.04  Apk-pk  to  52.14  Apk-Pk  following  the  loss  of  the  AC  grid  source,  which  the  SC  and 
associated  circuitry  are  well-suited  to  handle.  This  represents  an  additional  2.49  Apk-pk 
increase  in  SC  current  ripple  when  compared  to  the  “BPF  only”  scenario  but  allows  the 
removal  of  an  additional  10.06%  of  the  120-Hz  battery  bank  ripple  current,  bringing  the 
total  removal  to  99.99%  when  compared  to  the  “SC  Disconnected”  scenario.  A  close-up 
view  of  the  final  0.001  Apk-pk  battery  current  ripple  is  displayed  in  Figure  34. 
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Figure  34.  Battery  Current  120-Hz  Ripple,  Loss  of  AC  Grid  Source- 
SC  Bank  Connected  with  BPF  and  LPF 
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IV.  LABORATORY  IMPLEMENTATION  AND 
EXPERIMENTAL  MEASUREMENTS 


The  focus  of  the  laboratory  experimentation  for  this  thesis  was  to  show  that  the 
SC  bank  of  a  HESS  supplying  a  microgrid  can  be  used  to  remove  the  120-Hz  ripple  from 
the  battery  current,  thereby  extending  battery  life.  It  was  decided  that  certain  changes 
could  be  made  to  the  system  implementation  used  in  simulation  which  would  allow  the 
system  to  be  built  using  a  single  XILINX  Field  Programmable  Gate  Array  (FPGA)  while 
still  providing  confirmation  120-Hz  current  ripple  removal  from  the  battery  bank.  While 
the  laboratory  implementation  described  below  is  not  an  exact  physical  replication  of 
simulations  performed  for  this  thesis,  it  is  sufficient  to  show  that  SCs  can  be  used  to 
remove  the  120-Hz  current  ripple  from  the  battery  bank.  This  experimentation 
methodology  aims  to  mitigate  the  risk  (cost  and  resources)  of  building  a  complete 
physical  implementation  of  the  simulated  circuit  and  performing  rigorous  testing  by 
performing  limited  testing  that  still  validates  system  functionality. 

A.  DIFFERENCES  BETWEEN  SIMULATION  AND  PHYSICAL  SYSTEM 

The  differences  between  the  testing  performed  during  simulation  and  the  testing 
performed  with  the  physical  circuitry  built  in  the  laboratory  are  highlighted  in  this 
section.  In  each  subsection  below,  a  difference  between  the  two  test  methodologies  is 
identified  and  an  explanation  for  why  this  difference  does  not  affect  the  validity  of 
laboratory  analysis  is  provided. 

1.  Use  of  BPF  Only 

In  simulations,  a  BPF  was  used  to  isolate  the  120-Hz  DC  link  ripple  and  use  this 
ripple  to  increase  the  SC  boost  current,  removing  approximately  90%  of  the  120-Hz 
current  ripple  from  the  battery  bank.  A  FPF  was  added  to  the  boost  controller  and  used  to 
remove  additional  battery  current  ripple,  bringing  the  total  ripple  current  removal  to 
99.99%  with  the  use  of  both  filters.  While  this  BPF  with  FPF  configuration  produced  the 
largest  reduction  in  battery  ripple  current,  the  parallel  FPF  and  BPF  circuitry  (with 
associated  feedback)  is  more  difficult  to  build  and  program  in  the  laboratory  than  a  BPF 
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isolating  only  the  120-Hz  ripple  and  using  it  to  create  the  SC  boost  current.  By  using  a 
BPF  only,  the  reduction  in  battery  ripple  current  that  is  possible  with  the  addition  of  an 
SC  bank  can  still  be  validated  without  the  added  difficulty  of  programming  the  additional 
LPF,  so  a  LPF  was  not  used  for  laboratory  experimentation. 

2.  BPF  Gain  Reduced 

The  BPF  gain  K  used  for  simulations  was  set  at  30,  which  produced  SC  bank 
currents  of  up  to  30.0  A  in  the  laboratory.  This  current  value  is  beyond  the  steady-state 
current-carrying  capabilities  of  the  FPGA  and  associated  circuitry.  As  a  result,  the  BPF 
gain  was  reduced  to  K  =  2  so  that  useable  results  could  be  obtained  without  the 
possibility  of  overloading  laboratory  equipment. 

3.  SC  Voltage  Variation 

For  the  purpose  of  simulation,  it  was  necessary  to  assign  the  SC  bank  a  starting 
voltage  at  the  beginning  of  each  simulation.  While  the  value  of  this  voltage  affects  the 
battery  ripple  current,  it  was  decided  that  battery  ripple  current  reduction  could  be 
achieved  and  demonstrated  as  long  as  each  simulation  used  the  same  initial  SC  bank 
voltage.  An  arbitrary  starting  value  of  56.0  V  (the  steady-state  rated  voltage  limit  of  the 
130.0  F  Maxwell  SC  bank  that  was  used  in  the  laboratory)  was  chosen.  Additionally, 
although  the  SC  bank  voltage  is  changing  constantly  during  laboratory  experimentation 
based  on  DC  bus  voltage  and  SC  bank  current,  this  voltage  was  held  constant  at  56.0  V 
during  all  simulations.  This  made  it  easier  to  focus  the  output  plots  on  the  effect  of  using 
the  boost  controller  to  control  current  drawn  from  the  SC  bank  to  reduce  battery  ripple 
current  and  was  justified  in  that  SC  hank  voltage  is  not  expected  to  change  significantly 
over  the  course  of  a  short  simulation  (on  the  order  of  seconds).  In  laboratory 
experimentation,  the  SC  bank  voltage  changed  slightly  over  the  course  of  each 
experiment,  but  the  ability  of  the  SC  bank  to  provide  the  120-Hz  ripple  current  was  still 
clearly  visible. 
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4.  Supercapacitor  Current  Monitored  instead  of  Battery  Current 

Due  to  the  EMS  configuration  in  the  laboratory  and  the  space  available  for  probe 
placement,  it  was  more  convenient  to  measure  and  display  iscboost  than  i boost-  In  order  to 
provide  a  means  of  verifying  that  the  120-Hz  current  ripple  was  being  removed  from  the 
battery  bank  by  the  addition  of  the  SC  bank  and  BPF,  SC  hank  current  was  plotted  with 
and  without  the  use  of  a  BPF.  When  using  the  BPF,  the  120-Hz  ripple  that  was  removed 
from  battery  bank  current  was  directly  transferred  to  the  SC  bank,  so  demonstrating  the 
appearance  of  this  ripple  in  the  SC  bank  current  was  sufficient  to  prove  that  the  current 
ripple  was  removed  from  the  battery  bank. 


5.  Minimum  Usable  SC  Bank  Voltage 

For  all  simulations,  the  commanded  iscboost  and  iboost  currents  derived  by  the  boost 
controller  were  directly  applied  to  the  system  without  the  simulation  of  BBC  IGBT 
gating.  This  configuration  ignored  an  important  requirement  placed  on  the  physical 
system,  as  shown  in  [27],  which  is  that  the  SC  bank  BBC  must  be  maintained  in 
Continuous  Conduction  Mode  (CCM)  in  order  to  preserve  the  expected  converter  voltage 
transfer  function 


(19) 


where  Vc  is  the  voltage  across  the  BBC  output  capacitor  ( y/,„v  capacitance),  E  is  the  BBC 
input  voltage  (SC  bank  voltage),  and  D  is  the  IGBT  duty  cycle,  or  the  percentage  of  each 
switching  cycle  during  which  the  IGBT  is  off.  The  FPGA  circuitry  was  programmed 
assuming  that  this  transfer  function  was  valid  and  is  unable  to  properly  maintain  and 
control  system  parameters  (such  as  iscboost)  if  the  transfer  function  were  invalid.  As  the  SC 
bank  discharges  to  lower  voltages,  the  BBC  IGBT  duty  cycle  increases  to  maintain  iscboost 
equal  to  the  boost  controller  commanded  value  until  the  BBC  enters  Discontinuous 
Conduction  Mode  (DCM)  at  some  minimum  usable  SC  bank  voltage  value.  Once  this 
minimum  voltage  value  is  reached,  IGBT  duty  cycle  increases  are  no  longer  able  to 
maintain  the  commanded  iscboost  current.  In  practice,  SC  bank  voltage  dropped  as  low  as 
approximately  5.0  V  before  the  BBC  entered  DCM  and  data  became  invalid. 
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B.  LABORATORY  SETUP  AND  EQUIPMENT  USED 

The  EMS  was  built  in  the  laboratory  and  is  depicted  in  Figure  35  with  some  of  the 
major  components  labeled. 


Figure  35.  Laboratory  Setup  with  Major  Components  Labeled 

In  addition  to  the  major  components  depicted  in  Figure  35,  the  experiments  also 
required  a  laptop  with  XILINX  software  installed  to  program  and  interface  with  the 
FPGA  and  a  200  V  AC  power  supply  (for  grid  simulation). 

C.  DIGITAL  BPF  DESIGN 

The  laboratory  experimentation  involved  programming  a  XILINX  FPGA  to 
perform  all  of  the  functions  of  the  EMS  SIMULINK  model.  This  code  was  written  and 
provided  by  Naval  Postgraduate  School  professor  Alexander  L.  Julian,  Ph.D.  In 
simulations,  the  continuous-time  Laplace  domain  transfer  function  for  the  120-Hz  BPF, 
given  in  (18),  was  used.  The  first  task  in  designing  the  equivalent  digital  120-Hz  BPF  for 
programming  into  the  FPGA  was  to  convert  (18)  into  a  discrete-time  z-domain  difference 
equation.  In  the  MATLAB  script  provided  in  Appendix  A,  the  bilinear  transform  function 
of  MATLAB  was  first  used  to  find  the  coefficients  of  the  second-order  BPF  z-domain 
transfer  function,  producing 
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(20) 


F(z)  _  z2  (0. 1 125)  +  z(O.OOOO)  +  (-0. 1 125) 
X(z)  ~  z2(l  .0000)  +  z(-l  .9868)  +  (0.9925)  ’ 

where  Y(z)  is  the  filter  output,  X(z)  is  the  filter  input,  and 


z  = 


(21) 


where  /  is  operating  frequency  and  fs  is  sampling  frequency.  A  sampling  frequency  of 
10.0  kHz  was  chosen.  Rearranging  (20)  into  the  z-domain,  BPF  difference  equation 
produces 

Y(z)  =  (0.11 25) X(z)  +  z“2(-0. 1 125)X(z)  -  z’1  (- 1 ,9868)F(z)  -  z'2 (0.9925)F(z) .  (22) 

To  verify  that  (22)  defines  the  operation  of  a  120-Hz  digital  BPF  with  the  desired 
frequency  response,  the  MATLAB  script  in  Appendix  A  was  plotted.  From  Figure  36,  it 
can  be  seen  that  120-Hz  is  in  the  center  of  filter  pass  band. 


Figure  36.  Digital  120-Hz  BPF  Response 

By  inspecting  (22),  we  can  see  that  the  equivalent  digital  circuit  for  the  desired  BPF  can 
be  built  with  digital  multipliers,  adders,  subtractors,  and  delay  modules.  The  equivalent 
digital  120-Hz  BPF  is  as  shown  in  Figure  37. 
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Figure  37.  Equivalent  Digital  120-Hz  BPF 


The  SIMULINK  model  in  Figure  37  was  used  to  verify  that  (22)  produced  a  true 
equivalent  circuit  to  the  desired  120-Hz  analog  BPF  by  passing  a  120-Hz  sine  wave  into 
both  the  analog  and  digital  fdter  and  plotting  the  output  to  produce  Figure  38.  The  digital 
filter  output  is  shifted  downward  by  0.001  in  Figure  38  so  that  both  waveforms  are 
visible,  but  the  output  waveforms  are  otherwise  exactly  the  same  (as  expected). 
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Figure  38.  BPF  Output  For  Sine  Wave  Input,  Analog  and 

Equivalent  Digital  Filter 
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D.  LABORATORY  TESTING  AND  RESULTS 


The  EMS  DC  bus  controller  was  programmed  to  maintain  the  SC  bank  current  at 
a  constant  3.0  A  (charging  or  discharging)  when  the  BPF  was  not  being  employed  to 
allow  the  removal  of  battery  bank  120-Hz  ripple  current.  SC  charging  and  discharging 
operations  were  conducted  both  with  and  without  the  BPF,  and  source  voltage  v*  was 
plotted  with  SC  current  iscboost  in  each  case.  All  experiments  were  performed  (and  plots 
captured)  with  the  SC  bank  voltage  at  14.0  V. 

1.  SC  Charging 

When  in  charging  mode,  the  SC  bank  BBC  operates  as  a  buck  converter  in  CCM. 
Current  flowing  into  the  SC  bank  during  a  charge  is  arbitrarily  defined  as  the  negative 
direction  of  current,  based  on  current  sensor  placement.  The  EMS  was  placed  in  charging 
mode,  increasing  SC  bank  voltage  from  an  initial  0.0  V  to  14.0  V,  at  which  point  v.s  and 
iscboost  were  recorded.  Without  the  BPF,  vs  and  iscboost  are  as  plotted  in  Figure  39.  As 
expected,  average  SC  bank  current  is  approximately  -3.0  A  (with  a  slight  offset  due  to  a 
small  amount  of  sensor  inaccuracy),  and  there  is  no  120-Hz  ripple  present  in  iscboost- 
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Figure  39.  SC  Charging,  No  BPF 
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A  close-up  view  of  the  SC  bank  current  in  Figure  39  is  provided  in  Figure  40.  In 
this  figure,  it  is  clear  that  the  SC  bank  current  ripple  in  Figure  39  is  caused  by  the 
combination  of  the  15-kHz  IGBT  switching  frequency  and  duty  cycle. 
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Figure  40.  SC  Current  Ripple  Caused  By  Duty  Cycle  and 
15  kHz  Switching  Frequency 


Connecting  the  BPF  to  the  boost  controller,  the  DC  bus  controller  adjusts  SC  bank 
IGBT  duty  cycle  to  provide  the  iscboost  current  necessary  to  remove  the  120-Hz  ripple 
current  from  the  battery  bank.  In  this  configuration,  plots  of  vs  and  iscboost,  shown  in 
Figure  41,  clearly  demonstrate  that  the  120-Hz  ripple  was  removed  from  the  battery  bank 
and  placed  on  the  SC  bank.  This  ripple  has  an  amplitude  of  approximately  3.1  Apk-pk- 
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Figure  41.  SC  Charging,  with  BPF 
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2. 


SC  Discharging 


When  in  discharging  mode,  the  SC  bank  BBC  operates  as  a  boost  converter. 
Without  the  BPF,  v5  and  iscboost  are  as  plotted  in  Figure  42.  As  expected,  average  SC  bank 
current  is  approximately  3.0  A,  and  there  is  no  120-Hz  ripple  present  in  the  SC  hank 
current. 
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Figure  42.  SC  Discharging,  No  BPF 

Similar  to  the  results  obtained  during  the  “SC  Charging”  experiment,  connecting 
the  BPF  to  the  boost  controller  causes  the  DC  bus  controller  to  adjust  SC  bank  IGBT  duty 
cycle  to  provide  the  iscboost  current  necessary  to  remove  the  120-Hz  ripple  current  from  the 
battery  bank.  In  this  configuration,  plots  of  vs  and  iscboost,  shown  in  Figure  43,  again 
demonstrate  that  the  120-Hz  ripple  was  removed  from  the  battery  bank  and  placed  on  the 
SC  bank.  This  ripple  has  an  amplitude  of  approximately  2.3  Apk-pk. 
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Figure  43.  SC  Discharging,  with  BPF 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 
FOR  FUTURE  WORK 


The  work  conducted  throughout  this  thesis  provides  many  interesting  conclusions, 
as  well  as  room  for  potential  follow-on  research  and  experimentation.  In  this  section, 
some  of  the  more  pertinent  results  are  detailed  in  the  sincere  hope  that  future  students  and 
researchers  will  pick  up  the  mantle  and  continue  to  develop  this  system  into  a  useable 
end-product,  delivering  concrete  benefits  for  our  Armed  Forces  and  the  power  industry. 

A.  CONCLUSIONS 

The  primary  research  objective  of  this  thesis  was  to  show  through  both  simulation 
and  laboratory  experimentation  that  when  using  an  EMS  on  a  microgrid  with  a  single¬ 
phase  VSI,  the  battery  bank  ripple  current  caused  by  the  second-order  harmonic  DC  link 
voltage  ripple  can  be  filtered  out  by  adding  an  SC  bank  to  the  system  with  a  BPF  and 
boost  controller  circuitry.  In  the  accomplishment  of  this  objective,  several  specific 
conclusions  were  uncovered. 

1.  The  Benefits  of  Employing  an  EMS  Are  Clear 

The  simulations  and  laboratory  experimentation  performed  for  this  thesis  have 
clearly  shown  that  an  EMS  can  be  programmed  to  control  the  operation  and  utilization  of 
DR,  allowing  different  power  sources  and  energy  storage  methods  to  provide  services  to 
a  microgrid  in  ways  that  are  not  possible  without  an  EMS.  Specifically,  the 
characteristics  of  power  flow  to  and  from  different  HESS  energy  storage  devices  can  be 
individually  controlled  so  that  each  device  operates  more  efficiently.  This  is  especially 
important  when  dealing  with  batteries  as  an  energy  storage  method,  since  unnecessary 
cycling  caused  by  harmonic  ripples  in  DC  bus  voltage  can  quickly  degrade  battery 
capacity  and  lead  to  increased  energy  costs  when  battery  hanks  require  replacement. 
Additionally,  simulations  confirmed  that  the  microgrid  is  able  to  increase  power  system 
redundancy  by  maintaining  a  stable  power  supply  to  system  loads  when  the  grid  source  is 
lost,  which  is  highly  desirable  for  systems  powering  vital  loads.  The  optimization  of 
microgrid  parameters  is  achieved  without  any  negative  effect  or  need  for  control  of  the 
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grid  source,  making  the  microgrid  an  excellent  control  method  for  portable  systems 
requiring  a  stable,  redundant  power  source. 

2.  SCs  Can  Be  Used  to  Remove  Battery  Bank  120-Hz  Current  Ripple 

The  ability  to  remove  the  battery  bank  120-Hz  current  ripple  by  adding  an  SC 
bank  with  a  BPF  and  boost  controller  was  clearly  demonstrated  in  this  thesis.  In  two 
separate  simulation  scenarios,  this  additional  system  was  shown  to  reduce  battery  hank 
ripple  current  by  almost  90%  when  compared  to  operations  without  the  SC  bank  and 
BPF.  In  laboratory  experimentation,  the  SC  bank  current  was  measured  to  confirm  that 
the  120-Hz  current  ripple  was  removed  from  the  battery  bank  as  desired.  Even  at  the 
relatively  low  loading  and  current  values  achievable  in  the  laboratory  (based  on  system 
component  limitations)  the  removal  of  more  than  3.0  Apk_pk  of  ripple  current  from  the 
battery  bank  was  demonstrated.  While  not  explored  in  this  thesis  due  to  time  constraints, 
it  can  be  inferred  that  this  control  methodology  can  be  extended  to  other  harmonics  of 
concern  in  microgrids  employing  slightly  different  architectures  (such  as  three-phase  AC 
source  systems). 

3.  Almost  100%  of  Battery  Bank  120-Hz  Current  Ripple  Can  Be 
Removed  with  a  BPF  and  a  LPF 

As  a  secondary  research  objective,  the  possibility  of  further  reducing  battery  hank 
ripple  current  by  adding  an  LPF  to  the  boost  controller  was  explored  in  this  thesis.  Two 
different  simulations  provided  confirmation  that  a  battery  ripple  current  reduction  of  up 
to  99.99%  is  achievable  when  the  LPF  is  added  to  the  EMS  and  properly  tuned  as 
compared  to  a  system  employing  a  battery  bank  as  the  only  energy  storage  method.  This 
relatively  simple  addition  to  the  boost  controller  functionality  enables  the  system 
designers  to  achieve  an  almost  complete  removal  of  second-order  harmonic  ripple  current 
from  the  battery  bank,  bringing  the  system  close  to  maximum  efficiency. 

B.  RECOMMENDATIONS  FOR  FUTURE  WORK 

As  with  any  research  endeavor,  while  working  toward  the  research  objectives  for 
this  thesis,  the  topics  explored  raised  questions  and  brought  to  light  opportunities  for 
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future  work.  One  of  the  many  benefits  of  working  with  the  EMS  in  simulations  and  in  the 
laboratory  is  that  there  is  both  an  endless  universe  of  possible  expansion  areas  for  this 
research  and  a  great  need  for  this  expansion.  As  stated  in  the  introduction,  energy  costs 
are  an  enormous  burden  throughout  the  DOD;  anything  that  can  be  done  to  reduce  these 
costs  will  have  a  direct  positive  impact  on  our  nation.  The  recommendations  for  future 
work  focus  on  moving  closer  to  translating  the  current  system  into  a  useable  end-product 

1.  Test  with  a  Higher  Voltage  Capacity  SC  Bank  and  More  Robust  EMS 

The  laboratory  experimentation  for  this  thesis  was  limited  by  the  size  of  the  SC 
bank  and  the  voltage  and  current  limitations  of  EMS  circuitry.  For  example,  it  was  not 
possible  to  increase  the  battery  current  to  values  large  enough  to  allow  the  observation  of 
the  120-Hz  ripple  current,  although  it  was  clearly  shown  to  exist.  Using  an  SC  bank  with 
a  higher  voltage  capacity  allows  for  more  robust  experimentation,  such  as  the  physical 
implementation  of  the  two  simulations  performed  in  this  thesis  and  longer  experiments 
(due  to  the  larger  SC  energy  storage  capacity). 

2.  Add  Other  Power  Sources,  such  as  PV  Cells 

Adding  PV  cells  to  the  proposed  EMS  configuration  would  provide  yet  another 
means  by  which  system  efficiency  can  be  optimized.  With  the  increasing  prevalence  of 
solar  systems  used  as  alternative  sources  of  energy,  it  would  be  prudent  to  explore  how  a 
PV  cell  could  best  be  integrated  into  this  system. 

3.  Optimize  SC  Bank  Efficiency 

While  SCs  are  capable  of  operating  through  many  more  cycles  than  a  typical 
battery  before  requiring  replacement,  their  life  is  also  decreased  when  the  battery  bank 
current  ripple  is  shifted  to  the  SC  bank.  In  [28],  Abeywardana  et  al.  explore  a  current- 
feedback  method  that  reduces  SC  bank  harmonic  ripple  current  to  optimize  SC  bank  life. 
The  work  performed  in  this  thesis  can  be  expanded  to  include  this  method  or  another 
means  of  reducing  SC  bank  current  ripple  to  maximize  overall  system  efficiency  and 
cost-effectiveness. 
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APPENDIX  A.  MATLAB  SCRIPTS 


%%  LCDR  Gabriel  Hernandez,  Thesis 
%  With  assistance  from: 

%  NPS  Prof.  Giovanna  Oriti 
%  NPS  Prof.  Alexander  Julian 
%  NPS  Prof.  Roberto  Cristi 

'O 

%  Create  input  values  for  SIMULINK  model 

clear  all; 
close  all; 

%%  Establish  variable  values 
Kp_v  =  .00005; 

Ki_v  =  .05; 

sw_freq  =  15000; 

vo_ref  =  120*sqrt (2) *2/pi; 

vbus_ref  =  200;  %  Reference  DC  bus  voltage 

Lin  =  400e~6;  %  Includes  the  leakage  inductance  of  the  60Hz  xfmr, 

which  is  3mH 

Lfilter  =  (470+470+220) *le-6;  %  Includes  the  leakage  inductance  of  the 

60Hz  xfmr,  which  is  3mH 
Cin  =  12e-6; 

PWM_mode  =  1 ; 

Vbattery  =  100; 

Lbattery  =  100e-6; 

Rbattery  =  .1; 

Cbattery  =  1000e-6; 

Cbus  =  1000e-6; 

Lpv  =  5000e-6; 

Cpv  =  4000e-6; 
delta_I  =  .1; 
deltav  =  .004975; 

Isc  =  4; 

Rp  =  1; 

Rs  =  .05; 

Io  =  le-10; 

%%  PV  cell  calculations 

Vd  =  zeros ( 1 , 122 )  ; 

IPV  =  zeros ( 1 , 122 ) ; 

VPV  =  zeros  (1, 122)  ; 

Vd ( 1 )  =  .02; 

IPV(l)  =  Isc-Io* (exp (38 . 9*Vd(l) ) -1) -Vd(l) /Rp; 

VPV(l)  =  Vd(l) -IPV(l) *deltav; 
for  ii  =  2:122 

Vd(ii)  =  Vd (ii-1 ) tdeltav; 

IPV(ii)  =  Isc-Io* (exp (38 . 9*Vd (ii) ) -1) -Vd (ii) /Rp; 

VPV(ii)  =  Vd ( ii ) -IPV ( ii ) *deltav; 

end 

for  ii  =  1:122 
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IPV_rev ( ii )  =  IPV(123-ii); 

VPV_rev(ii)  =  VPV(123-ii); 

end 

IPV_LU  =  [1 : 122] /122*4; 

VPV_LU  =  [1 : 122] *0; 

IPV_LUindex  =  [1:122]; 
nn  =  2 ; 

for  ii  =  1:121 

while  IPV_rev (nn) <IPV_LU (ii) 
nn  =  nn+1; 

end 

slope  =  (VPV_rev (nn) -VPV_rev (nn-1 ) ) / ( IPV_rev (nn) -IPV_rev (nn-1 ) )  ; 

VP V_LU ( i i )  =  slope* (IPV_LU (ii) -IPV_rev (nn-1) ) +VPV_rev (nn-1)  ; 

end 

Power_PV  =  VPV . * IPV ; 

%%  Values  most  often  adjusted 
%  Simulation  start/stop  and  step  times 

pwm_on  =0;  %H  bridge  is  PWM  modulated  when  this  variable  is  1 
tstop  =  17; 

if  pwm_on  ==  0 

tstep  =  15e-6; 

else 

tstep  =  10e-6; 

end 

%  Turn  source  voltage  on/off 

source_on  =20;  %  time  when  source  connects 

source_off  =  25;  %  time  when  source  disconnects 

%  Added  portion  for  supercapacitor 

vsc  =  56;  %  Voltage  limit  of  supercapacitor  is  56V.  In  laboratory,  will 

start  from  complete  discharge 

rsc  =  8.1e-3; 

lsc  =  1000e-6; 

esc  =  130; 

sc_connection_status  =  1;  %"1"  =  SC  bank  connected,  "0"  =  SC  bank 
disconnected.  BPF  is  connected  with  SC  Bank 

lpf_connection_status  =  1;  %"1"  =  LPF  connected,  "0"  =  LPF  disconnected 

%  BPF  for  120Hz  oscillation 
K  =  30; 

bpf_freq  =  120; 
bpf_wo  =  2*pi*bpf_f req; 
lpf_freq  =  .1; 
lpf_wc  =  2*pi*lpf_f req; 

Q  =  10; 

%  Allow  for  a  load  change  during  operation 
Rload  =  200; 

Rload2  =  10; 
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Rload_applied  =  0; 

Rload2_applied  =  5; 

%  Change  EMS  reference  current 
iref_change_at_time  =  20; 
iref_initial_value  =  1; 
iref_f inal_value  =  4; 

%%  BPF  Design 

%  Calculate  BPF  Z-transform  coefficients  using  bilinear  transformation 
%  Gives  the  Z-domain  transfer  function  required  for  the  equivalent 
%  digital  filter  to  the  analog  filter  governed  by  the  chosen  S-domain 
%  transfer  function 

fs  =  10e3;  %  Assume  10kHz  sampling  frequency 

[numd,dend]  =  bilinear ([ (K)  0 ]  ,  [  (Q/bpf_wo)  1  (Q*bpf_wo) ] , f s) ; 

%  Verification  of  BFP  frequency  response. 

low_freq  =  100;  %  Lowest  plot  frequency 
high_freq  =  140;  %  Highest  plot  frequency 
points  =  1000; 

spacing  =  (high_freq  -  low_freq) / (points  -  1);  %  Spacing  between 
frequencies 

for  i  =  1 :  points 

freq_val(i)  =  low_freq  +  (i-1) *spacing; 

Z  =  exp(j  *  2  *  pi  *  freq_val(i)  /  fs);  %  The  definition  of  Z 
H_z  =  ( ( Z^2 ) *numd ( 1 )  +  Z*numd(2)  +  numd(3))/... 

(  (ZA2) *dend(l)  +  Z*dend(2)  +  dend(3)); 
mag(i)  =  abs(H_z); 
phase  =  angle (H_z) *180/pi; 

end 

%  Plot  BPF  output 
figure  ( 12 ) 

plot (freq_val,  20*logl0 (mag) ) ; 

xlabel ( 'Freq  (Hz) ' ) ; 

ylabel (' Filter  response  (dB) ') 

set (gca, ' FontSize ' , 12 ) 

grid  on 

print (gcf, ' -dtif f ' -r350 ' , 'BPF  output ' ) 
savefig('BPF  Output') 
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%%  LCDR  Gabriel  Hernandez,  Thesis 
%  With  assistance  from: 

%  NPS  Prof.  Giovanna  Oriti 
%  NPS  Prof.  Alexander  Julian 
%  NPS  Prof.  Roberto  Cristi 

'o 

%  Plot  SIMULINK  waveforms 

close  all 

v_ems=data_out ( : , 1 ) ; 
i_ems=data_out ( : ,  3) ; 
i_s=data_out ( :  ,  4 )  ; 

i_ref =data_out2 ( :  ,  1)  ; 

v_ref =data_out2 ( :  ,  2 )  ;  %  Reference  voltage  when  grid-connected 

%%  data_battery  -  Battery  MUX  output 
v_bus=data_battery ( :  ,  1)  ; 
v_batterycap=data_battery ( : ,  2 )  ; 
i_boost=data_battery ( : ,  3) ; 

i_batt=data_battery ( : ,  4 ) ;  %  actual  battery  current 

i_bus=data_battery ( : ,  5)  ; 

i_emsdc=data_battery ( : , 6 ) ;  %  into  DC  bus... called  "ibus"  in  simulink 

model 

%%  data_sc  -  SC  MUX  output 

v_sc=data_sc ( : , 1 ) ;  %  Voltage  across  sc 

i_scboost=data_sc ( : , 2 ) ;  %  Current  in  scbus 

i_scbus=data_sc ( : , 3 ) ;  %  supercapacitor  boost  current 

%%  data_out3  -  MUX  output  from  physical  system 
real_and_ind_iload=data_out3 ( : , 1) ; 
iload= (data_out3 ( : ,  2 ) )  ; 
v_ac  =  (data_out3 ( : ,  3) )  ; 

SC_power_out  =  v_sc . *i_scboost; 
batt_power_out  =  v_batterycap . *i_boost; 

t_ss=round ( . 1/tstep)  ; 
one_cycle=round (1/60/tstep)  ; 

vac_ss=v_ems (t_ss-2*one_cycle : t_ss-one_cycle ) ; 
vac_spect=f f t (vac_ss) / (one_cycle/2 ) ; 
f req_vac= [ 0 : one_cycle ] *  60 ; 
v_inv=data_out ( :  ,  2 )  ; 

%%  data_out4  -  Boost  controller  transformation  verification  output 

sine_in  =  data_out4 ( : , 1 ) ; 
analog_bpf_op  =  data_out4 ( : , 2 ) ; 
discrete_bpf_f unc_op  =  data_out4 ( : , 3) ; 
h_domain_built_bpf_op  =  data_out4 ( : , 4) ; 

%  Verification  plot  for  boost  controller 
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figure (11) 

plot (bpf_time, analog_bpf_op, ' r-. ') 
hold  on 

xlabel('Time  (s) ') 

ylabel('BPF  output  for  sine  wave  input') 
plot (bpf_time, h_domain_built_bpf_op- . 001 , ' b ' ) 
legll  =  ... 

legend (' Output  from  analog  f ilter ',' Output  from  equivalent  digital 
f ilter- . 00 1 ' , 'location ' , ' best ' ) ; 
figll_leg  =  f indob j ( legl 1 ,' type ',' text ')  ; 
set (f igll_leg,  ' FontSize ' , 12 ) 
set (gca, ' FontSize ', 12 ) 
grid  on 
hold  off 

print (gcf,  ' -dtiff ',' -r350 ',' Boost  Controller  Verification') 


%%  Plot  figures 


figure ( 1 ) ; 


plot (time, i_scboost,  ' g ' ) 
hold  on 

plot (time, i_batt, ' r ' ) 
xlabel('Time  (s) ') 
ylabel (' Current  (A) ') 

legl  =  legend (' i_{ scboost }',' i_{batt }',' location ',  'best ') ; 

figl_leg  =  findobj (legl,  'type ',  'text ') ; 

set (f igl_leg, ' FontSize ', 12 ) 

set (gca, ' FontSize ', 12 ) 

grid  on 

print (gcf ,' -dtiff ',' -r350 ',' isc_boost  and  i_batt ' ) 
savef ig ( ' isc_boost  and  i_batt ' ) 


figure  (2 ) 


plot (time, v_bus, ' b ' ) ; 
xlabel('Time  (s) ') 
ylabel (' v_{bus }  (V) ') 

leg2  =  legend (' v_{bus }',' location ',' best ') ; 
fig2_leg  =  f indob j (leg2 ,' type ',' text ') ; 
set ( f ig2_leg,  ' FontSize ', 12 ) 
set (gca, ' FontSize ', 12 ) 
grid  on 

print (gcf , ' -dtiff ' , ' -r350 ' , ' v_bus ' ) 
savef ig ( ' v_bus ' ) 
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figure ( 3 ) 


plot (time, i_batt,  ' r ' ) 
xlabel('Time  (s) ') 
ylabel (' Current  (A) ') 

leg3  =  legend (' i_{batt location best ')  ; 
fig3_leg  =  f indob j ( legl ,  ' type text ') ; 
set ( f ig3_leg, '  FontSize ', 12 ) 
set (gca, ' FontSize 12 ) 
grid  on 

print (gcf ,  ' -dtif f ' ,  ' -r350 '  ,  ' i_batt ' ) 
savef ig ( ' i_batt ' ) 

%%  FFT  plot  of  i_ems  to  show  60  and  180  Hz  components 


figure  (10) 


one_period=round (1/60/tstep) ; 

v_ems_ss=v_ems (2/tstep : 2/tstep+one_period, 1 ) ;  %  data  at  t=4  seconds 
v_ems_spect=f f t (v_ems_ss) / (one_period/2 )  ; 

plot (0, 20*logl0 (abs (v_ems_spect (1) ) /2) , 'bo ' ) ;  %  DC  component 

hold  on 

plot ( ( [2:40] -1) , 20*logl0 (abs (v_ems_spect (2:40))), ' ro ' ) ; 
hold  off 

xlabel (' Harmonic  Number'); 
ylabel ( ' i_{ems }  [dB]  ' )  ; 


60 


%%  LCDR  Gabriel  Hernandez,  Thesis 
"6 

%  With  assistance  from: 

%  NPS  Prof.  Giovanna  Oriti 
%  NPS  Prof.  Alexander  Julian 
%  NPS  Prof.  Roberto  Cristi 

'o 

%  Script  to  plot  laboratory  results 

%%  Clear  variables  and  command  window,  and  close  all  plots 
clc; 

clear  all; 
close  all; 

%%  Import  Excel  spreadsheet  data 

%  Charging,  with  120  Hz  ripple 

vsx_chl20  =  xlsread ( ' Tek _ CHI _ Wfm _ vs . csv ' ) ; 

iSCx_chl20  =  xlsread (' Tek_CH3_Wfm_iSC . csv ') ; 

%  Charging,  without  120  Hz  ripple 

vsx_chnol20  =  xlsread (' Tek_CHl_Wfm_vs_nol20 . csv ') ; 
iSCx_chnol20  =  xlsread (' Tek_CH3_Wfm_iSC_nol20 . csv ') ; 

%  Discharging,  with  120  Hz  ripple 

vsx_disl20  =  xlsread (' Tek_CHl_Wfm_vs_discharge . csv ') ; 
iSCx_disl20  =  xlsread (' Tek_CH3_Wfm_iSC_discharge . csv ') ; 

%  Discharging,  without  120  Hz  ripple 

vsx_disnol20  =  xlsread (' Tek_CHl_Wfm_vs_nol20_discharge . csv ') ; 
iSCx_disnol20  =  xlsread (' Tek_CH3_Wfm_iSC_nol20_discharge . csv ') ; 

%%  Plot  output  data 

%  Charging,  with  120  Hz  ripple 

len_chl20  =  length (vsx_chl20 ) ; 
time_vec_chl20  =  vsx_chl20 (15 : len_chl20, 1) ; 

time_vec_chl20  =  time_vec_chl20+ . 02 ;  %  Shift  the  time  axis  to  start  at 
t  =  0 

vs_chl20  =  vsx_chl20 (15 : len_chl20, 2) ; 
iSC_chl20  =  iSCx_chl20 (15 : len_chl20, 2) ; 
tmin_chl20  =  time_vec_chl20 ( 1 ) ; 
tmax_chl20  =  max (time_vec_chl20 ) ; 

figure  ( 1 ) ; 

plot (time_vec_chl20 , vs_chl20/ 20 , time_vec_chl20 , -iSC_chl20 ) ; 
xlabel('Time  (s) ') 

ylabel (' Volts  (V)  or  Current  (A)'); 
legl  =  ... 

legend ( ' v_s /2  0 ' ,  ' i_{  scboost }  ' ,  'Location ' ,  'Northeast ' ) ; 
figl_leg  =  f indob j ( legl ,' type ',' text ') ; 
set ( f igl_leg, ' FontSize ' , 12 ) 
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set (gca, ' FontSize ' , 12 ) 
grid  on; 

axis ( [tmin_chl20  tmax_chl20  -8  8]) 

print (gcf,  ' -dtiff -r350 Charging  with  120  Hz  ripple'); 

%  Charging,  without  120  Hz  ripple 

len_chnol20  =  length (vsx_chnol20 ) ; 
time_vec_chnol20  =  vsx_chnol20 (15 : len_chnol20, 1) ; 

time_vec_chnol20  =  time_vec_chnol20+ . 02 ;  %  Shift  the  time  axis  to  start 
at  t  =  0 

vs_chnol20  =  vsx_chnol20 (15 : len_chnol20, 2) ; 
iSC_chnol20  =  iSCx_chnol20 (15 : len_chnol20, 2) ; 
tmin_chnol20  =  time_vec_chnol20 (1)  ; 
tmax_chnol20  =  max (time_vec_chnol20 )  ; 

figure (2 ) ; 

plot (time_vec_chnol2  0 , vs_chnol2  0/ 20 ,  time_vec_chnol2  0 ,  -  iSC_chnol2  0 )  ; 
xlabel('Time  (s) ') 

ylabel (' Volts  (V)  or  Current  (A)'); 
leg2  =  ... 

legend ( ' v_s /2  0 ' ,  ' i_{  scboost }  ' ,  'Location ' ,  'Northeast ' ) ; 
fig2_leg  =  f indob j ( leg2 ,' type ',' text ') ; 
set ( f ig2_leg, ' FontSize ', 12 ) 
set (gca, ' FontSize ', 12 ) 
grid  on; 

axis ( [tmin_chnol20  tmax_chnol20  -8  8]) 

print (gcf ,' -dtiff ',' -r350 ',' Charging  without  120  Hz  ripple'); 

%  Discharging,  with  120  Hz  ripple 

len_disl20  =  length (vsx_disl20 ) ; 
time_vec_disl2  0  =  vsx_disl2  0 (15 : len_disl2  0,  1) ; 

time_vec_disl20  =  time_vec_disl20+ . 02 ;  %  Shift  the  time  axis  to  start 
at  t  =  0 

vs_disl20  =  vsx_disl20 (15 : len_disl20, 2) ; 
iSC_disl20  =  iSCx_disl20 (15 : len_disl20, 2) ; 
tmin_disl20  =  time_vec_disl20 ( 1 ) ; 
tmax_disl20  =  max (time_vec_disl20 ) ; 

figure ( 3 ) ; 

plot (time_vec_disl2  0 , vs_disl2  0/ 20 , time_vec_disl2  0 ,  -iSC_disl2  0 )  ; 
xlabel ( ' Time  (s) ' ) 

ylabel (' Volts  (V)  or  Current  (A)'); 
leg3  =  ... 

legend ( ' v_s /2  0 ' ,  ' i_{ scboost } ' ,  'Location ' ,  'Northeast '  ) ; 
fig3_leg  =  f indob j ( leg3, ' type ',' text ') ; 
set ( f ig3_leg,  ' FontSize ', 12 ) 
set (gca, ' FontSize ', 12 ) 
grid  on; 

axis ( [tmin_disl20  tmax_disl20  -8  8]) 

print (gcf ,  ' -dtiff ',' -r350 ',' Discharging  with  120  Hz  ripple'); 


%  Discharging,  without  120  Hz  ripple 
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Ien_disnol20  =  length (vsx_disnol20 ) ; 

time_vec_disnol20  =  vsx_disnol20 (15 : len_disnol20, 1) ; 

time_vec_disnol20  =  time_vec_disnol20+ . 02 ;  %  Shift  the  time  axis  to 
start  at  t  =  0 

vs_disnol20  =  vsx_disnol20 (15 : len_disnol20, 2) ; 
iSC_disnol20  =  iSCx_disnol20 (15 : len_disnol20, 2) ; 
tmin_disnol20  =  time_vec_disnol20 ( 1 )  ; 
tmax_disnol20  =  max (time_vec_disnol20 ) ; 

figure ( 4 )  ; 

plot (time_vec_disnol2  0 ,  vs_disnol2  0/ 20 ,  time_vec_disnol2  0 ,  -iSC_disnol2  0 )  ; 
xlabel('Time  (s) ') 

ylabel (' Volts  (V)  or  Current  (A)'); 
leg4  =  ... 

legend ( ' v_s /20 ' , ' i_{ scboost } ' , 'Location ' , 'Northeast ' ) ; 
fig4_leg  =  f indob j ( leg4 ,' type ',' text ') ; 
set ( f ig4_leg,  ' FontSize ' , 12 ) 
set (gca, ' FontSize ', 12 ) 
grid  on; 

axis ( [tmin_disnol20  tmax_disnol20  -8  8]) 

print (gcf,  ' -dtiff ',' -r350 ',' Discharging  without  120  Hz  ripple'); 

%  Additional  plot  to  show  effect  of  IGBT  gating 
figure ( 5 ) ; 

plot (time_vec_chnol20,  -iSC_chnol20,  '  r  '  )  ; 

xlabel('Time  (s) ') 

ylabel ( 'Current  (A) ' ) ; 

set (gca, ' FontSize ', 12 ) 

grid  on; 

axis ( [ . 0201  .0206  -4.5  -1]) 

print (gcf ,' -dtiff ',' -r350 ',' Charging  without  120  Hz  ripple'); 
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APPENDIX  B.  LABORATORY  OUTPUT  VALUES 


Simulation 

Condition 

SC  Bank  and 
Filter  Usage 

Battery 

Current 

Ripple 

Before 

Change 

(App) 

Average 

Battery 

Current 

Before 

Change 

(A) 

Grid 

connected, 

operator 

increases 

commanded 

current 

SC  bank 

disconnected 

1.006 

0.8264 

SC  bank 

connected 

with  120  Hz 
BPF  only 

0.1013 

0.82615 

SC  bank 

connected 

with  120  Hz 

BPF  and  LPF 

0.0001 

0.82615 

Loss  of  AC 
grid  source 
and 

SC  bank 

disconnected 

1.0063 

0.82645 

SC  bank 

connected 

with  120  Hz 

0.1013 

0.82615 

lr aribi Lion  lu 

islanding 

mode 

SC  bank 

connected 

with  120  Hz 

BPF  and  LPF 

IE- 04 

0.82395 

Islanding, 

step 

increase  in 

load 

SC  bank 

disconnected 

2.3783 

0.71485 

SC  bank 

connected 

with  120  Hz 
BPF  only 

0.2408 

0.714 

SC  bank 

connected 

with  120  Hz 

BPF  and  LPF 

0.0002 

0.714 

Battery 

Current 

Ripple 

After 

Change 

(App) 

Average 

Battery 

Current 

After 

Change 

(A) 

3.9485 

3.35845 

0.399 

3.4001 

0.0003 

3.35055 

17.0307 

14.60035 

1.7623 

14.81495 

0.0016 

14.8303 

17.0306 

14.6003 

1.7627 

14.81475 

0.0015 

14.82955 

89.93% 


89.93% 


89.88% 


99.99% 


99.99% 


89.89% 


89.65% 


89.65% 


99.99% 


99.99% 


SC  Bank 

SC  Bank 

Increase  in 

Current 

Current 

SC  Current 

Ripple 

Ripple 

Ripple 

Before 

After 

Caused  By 

Change 

Change 

Change 

(App) 

(App) 

(A) 

99.99% 


99.99% 


2.8753 

11.263 

8.388 

3.0483 

11.922 

8.874 

2.8767 

49.4784 

46.602 

3.0443 

52.1405 

49.096 

6.8321 

49.4786 

42.647 

7.2355 

52.1406 

44.905 

Specifics 


Change  at  10.0 
sec 

EMS  current 
changed  from 
1.0A  to  4.0A 
System  load  = 
lOohms 


Change  at  5.0 
sec 

Grid  source 
removed 
System  load  = 
lOohms 
EMS 

commanded 
current  =  1A 

Change  at  5.0 
sec 

System  load 
increased  from 
200  ohms  to  10 
ohms 
EMS 

commanded 
current  =  1A 
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